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Abstract 

The dynamics of electrons and atoms interacting with intense and 
ultrashort optical pulses presents an important problem in physics 
that cuts across different materials such as semiconductors and met- 
als. The currently available laser pulses, as short as 5 fs, provide a 
time resolution shorter than the dephasing and relaxation times in 
many materials. This allows for a systematic study of many-body 
effects using nonlinear optical spectroscopy. In this review article, we 
discuss the role of Coulomb correlations in the ultrafast dynamics of 
modulation-doped quantum wells and metal nanoparticles. We focus 
in particular on the manifestations of non-Markovian memory effects 
induced by strong electron-hole and electron-plasmon correlations. 
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1 Introduction and Scope of this Review 

One of the most exciting thrusts in Physics focusses on understanding the 
properties of systems possessing many degrees of freedom. The goal is to 
relate such properties to the interactions among the elementary excitations. 
Such many-body problems have been a driving force in many different sub- 
fields of Physics and Chemistry: condensed matter physics, nuclear physics, 
quantum chemistry, high energy physics, • • • In condensed matter physics, 
the study of many-body phenomena has led to many exciting discoveries: 
superconductivity, quantum Hall effect, Kondo effect, Mott transition, quan- 
tum phase transitions, Fermi Edge and X-ray Edge Singularities, ■ ■ ■. In 
most cases, such fascinating phenomena have been observed in equilibrium 
systems at low temperatures. In comparison, the role of many-body effects 
in the dynamics of non-equilibrium many-electron systems presents a less 
explored frontier. 
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Ultrafast nonlinear optical spectroscopy provides unique and powerful 
tools for studying the dynamics of many-body correlations. In such exper- 
iments, an intense optical pulse, the pump, with duration shorter than the 
scattering time or the period of the elementary excitations photoexcites the 
system from its ground state. The subsequent dynamics is then monitored 
by using a second pulse, the probe. By measuring the amplitude and the 
phase of the optical field emerging from the photoexcited sample, one can 
gain valuable insight into the interplay between many-body and quantum 
confinement effects during very short time scales. A good understanding of 
the role of such fundamental physics can lead to new ideas regarding the 
next generation of photonic and opto-electronic devices. These operate un- 
der intense excitation conditions in the technologically important Tbit/sec 
regime and are often based on semiconductor heterostructures and novel low- 
dimensional structures. Such considerations have spurred intense efforts in 
recent years aimed at understanding how to control the optical dynamics 
during sub-picosecond time scales (1 Tbit/sec = (1 ps)^^). In the highly 
non-equilibrium femtosecond regime, the dominant effects on the nonlinear 
optical dynamics come from the Pauli exclusion principle (Phase Space Fill- 
ing effects) and the Coulomb interactions among the optically-excited and 
Fermi sea carriers. The phonons also play an important role, especially dur- 
ing picosecond time scales. 

Following their photoexcitation, the electrons and holes undergo a num- 
ber of different relaxation stages as they scatter among themselves via the 
Coulomb interaction and with the phonons via the electron-phonon [e-p) in- 
teraction. Below we briefly discuss each of these relaxation stages. Initially 
(stage I, the coherent regime) the photoexcited electrons and holes have a 
well-defined phase relationship: since the photon momentum is negligible, 
they are restricted to have opposite momenta. Such a coherence is described 
by an optical polarization. In the very beginning, the photoexcited carri- 
ers interact via the bare Coulomb interaction. After time intervals of the 
order of the inverse plasmon frequency corresponding to the carrier density, 
the interactions become screened. The problem can then be approached in 
terms of interacting Fermi liquid quasi — particles [|^ . Even though in typical 
metals such screening occurs within an extremely short time interval ~ 1 
fs, in semiconductors the inverse plasma frequency is of the order of tens 
of femtoseconds and therefore the above regime is accessible experimentally 
by using ~ 5-10 fs optical pulses. If one tunes the pump frequency above 
the onset of linear absorption, real carriers are photoexcited, in which case 
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the duration of stage I is determined by the dephasing time, T2. In metals, 
the latter is typically of the order of a few femtoseconds, while in semicon- 
ductors it can be as long as a few picoseconds. The dephasing processes 
are suppressed if one tunes the pump frequency below the onset of linear 
absorption. Under such below-resonance excitation conditions, the coherent 
effects become dominant since only virtual carriers are photoexcited and thus 
dissipative scattering processes cannot occur. The time-dependence of the 
optical polarization is then determined by the pulse duration, which leads to 
extremely fast nonlinearities. 

Immediately after the e-/i phase coherence is lost, incoherent populations 
of electrons and holes are formed. These are initially described by non- 
thermal density distributions. Subsequently, these evolve into hot thermal 
distributions characterized by an electron temperature that can far exceed 
that of the lattice. Such a time evolution (stage II) as the carriers equi- 
librate among themselves is mainly induced by the e-e, e-h, h-h, and e-p 
interactions. The duration of this stage II is determined by the population 
relaxation time, Ti, which can be of the order of picoseconds. Subsequently, 
the e-p interaction leads to the equilibration of the electronic and lattice 
temperatures on a time scale of tens of picoseconds, followed by the return 
to thermodynamic equilibrium after hundreds of picoseconds via the recom- 
bination of the photoexcited carriers or via the transport of heat to the 
surroundings (stage III). 

A large body of experimental and theoretical work has focussed on dif- 
ferent aspects of the above relaxation. Such developments in undoped semi- 
conductors have been extensively described in a number of excellent review 
articles and books. Among the earlier reviews of interaction effects in highly- 
excited undoped semiconductors, we note those of Haug and Schmitt-Rink 
0], Stahl and Baslev [^, Zimmermann Schmitt-Rink et. al. and Haug 
and Koch . More recent reviews focussing on the femtosecond coherent dy- 
namics in undoped semiconductors include those by Mukamel [0 , Binder and 
Koch P), Shah 0, Haug and Jauho [jlO[, Axt and Mukamel [0], and Chemla 
P^ . We also note the collections of papers in Refs. [|T3], |T^, [T^. More 
recently, there has been a lot of interest in using ultrafast nonlinear spec- 
troscopy to study the dynamics of many-body effects in metals. The charge 
carrier dynamics and the coherent ultrafast phenomena in bulk metals, metal 
interfaces, and metal surfaces were recently reviewed by Petek and Ogawa 



The purpose of this review article is to provide a comprehensive picture 
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of recent experimental and theoretical advances in understanding femtosec- 
ond many-body phenomena in Fermi sea systems. As has been extensively 
discussed in the review articles and books listed above, even though a de- 



scription of the ultrafast dynamics based on few-level models ||T8|, [T^, |20 
has been very successful in atomic and molecular systems, it faces serious 
shortcomings when applied to semiconductors and metals. The development 
of microscopic non-equilibrium many-body theories has been found to be 
necessary for describing the femtosecond dynamics and has made this field 
an exciting area of research in condensed matter physics. In view of the 
large body of work, especially in undoped semiconductors, we will concen- 
trate here on specific aspects of the many-body physics that have not been 
extensively reviewed elsewhere. In particular, we will not touch on the im- 
pressive advances in undoped semiconductors and will attempt to provide 
a unified picture of the ultrafast dynamics in different confined Fermi sea 
systems. We will mainly focus on coherent effects due to the Coulomb in- 
teraction and will only briefly discuss the electron-phonon [e-p) interaction 
effects. The main body of this article discusses the recent breakthroughs 
in understanding the sub-picosecond dynamics in modulation-doped quan- 
tum wells (MDQW's) and metal nanoparticles, with particular emphasis on 
non-Markovian memory effects induced by carrier-carrier many-body corre- 
lations. Such dynamical effects are beyond the scope of the dephasing and 
relaxation time approximations. One should note here that a detailed quan- 
titative description of the experimental results presents a formidable task 
in the subpicosecond regime, which requires sophisticated numerical simula- 
tions as well as different approximations. This can make it difficult to discuss 
in a simple and intuitive way the important many-body physics conveyed by 
the experiment. Here we will focus on the simplest possible models that 
can capture the most important physics associated with the non-equilibrium 
Coulomb correlations and attempt to provide, to the degree possible, an 
intuitive explanation of the often delicate experimental features. 

Let us now briefly discuss how many-body correlations enter into the 
ultrafast nonlinear optical response. The most commonly used theoretical 
approach starts with the Heisenberg equations of motion of the one-particle 
density matrix, Pij{t), where the indices i,j = 1,2 label the conduction and 
valence band 0. The off-diagonal density matrix element, puit) = P(t), 
is the optical polarization, while the diagonal matrix elements, pii{t), corre- 
spond to the electron and hole occupation numbers (density distributions). 
Such equations of motion must be solved to obtain the polarization, which 
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in turn acts as the source term in the Maxwell equations that determine 
the signal measured in the experiment [0. The difficulty in treating the 
many-body interactions comes from the coupling of the one-particle den- 
sity matrix to many-particle correlation functions (higher density matrices). 
As is typical for many-body systems, this leads to an infinite hierarchy of 
coupled equations which cannot be solved exactly |lT|. One must there- 



fore devise a controlled truncation scheme. One possibility is to adopt the 
time-dependent Hartree-Fock approximation (HFA), which relies on the fac- 
torization of the many-particle correlation functions and only includes the 
two-particle interaction effects. Such an approach developed into the very 
successful Semiconductor Bloch Equations (SBE's) formalism (see e.g. Refs. 
2T|, P^, P^, El, p3, However, many recent experimental observations. 



reviewed e.g. in Ref. |12[, could not be explained within such a mean field 
approach and were attributed to four-particle and higher order Coulomb 
correlations. In order to develop a many-body theory that can capture such 
effects, it was noted by several different groups that, in undoped semicon- 
ductors, the above infinite density matrix hierarchy truncates if one adopts 
an expansion in terms of the optical fields 0, 0, If. This 

is due to the fact that, in the undoped semiconductor ground state, the 
conduction band is empty while the valence band is full. Such a method 
for treating the Coulomb correlations nonperturbatively is often referred to 
as the "dynamic controlled truncation scheme" . An alternative many-body 
approach has been based on diagrammatic expansions for Keldysh Green's 
functions, whose multiple-time dependence is typically eliminated by using 
the Kadanoff-Baym ansatz . 



Let us now discuss how non-Markovian memory effects enter into the 
density matrix equations of motion. Due to the many degrees of freedom 
that affect the dynamics, one usually distinguishes between a subsystem in- 
teracting directly with the optical fields (e.g. the photoexcited e-h pair) 
and a reservoir /bath consisting of all the other degrees of freedom (e.g. the 
phonons or the Fermi sea excitations). The polarization equation of motion 
can be cast in the form 

= l^Wlcoh + l^Wlscatt (1) 

where the ffist term describes the Hartree-Fock interactions among the co- 
herent e-h pairs, while the second term describes the many-body scattering 
processes among the carriers as well as between the carriers and the bath 
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excitations. The latter term can be expressed in the general form 

|^WIscatt = f_jmt,t')P{t') (2) 

where T{t,t') is a memory kernel (similar equations can also be written for 
the carrier distribution functions). The importance of the many-body cor- 
relations of interest here is that they determine the dependence of r(t, t') on 
the two times t and t' and therefore govern the memory effects. In systems 
where only times t' ^ t contribute to the above integral, one can apply the 
Markovian approximation, in which case 

|n^)lscatt-r(t)P(t). (3) 

In many cases, one can further approximate r{t) by a time-independent con- 
stant, the dephasing energy width ~ I/T2. For example, such an approxima- 
tion allows one to include, in a semi-phenomenological way, the effect of the 
scattering processes in the SBE's. However, when such a Markovian approx- 
imation breaks down, the many-body correlations can affect the qualitative 
ultrafast dynamics by inducing strong memory effects and non-exponential 
polarization decay. Such effects in undoped semiconductors are known to 
arise e.g. from exciton-exciton correlations or from e-p and carrier-carrier 
interactions in the quantum kinetic regime and have been extensively re- 
viewed e.g. in Refs. ||T2|, |10|]. Here we focuss on analogous effects whose 
origin lies in the non-perturbative Coulomb correlations between the pho- 
toexcited carriers and the cold electron Fermi sea in MDQW's and small 
metal nanoparticles. 

This review is organized as follows. First we present an overview of the 
experimental results on the ultrafast nonlinear optical dynamics in MDQW's 
(Section |^) and metal nanoparticles (Section H). In Section § we outline the 
formalism necessary for describing the coherent ultrafast nonlinear optical re- 
sponse of Fermi sea systems. In Section |^, the ultrafast nonlinear dynamics 
of the Fermi Edge Singularity (FES) in the coherent regime is discussed. In 
Section ^ the size-dependent dynamical screening of the Coulomb interaction 
in metal nanoparticles is discussed and the results for the quasiparticle scat- 
tering rates are reviewed. In Section ^ these results are used to describe the 
ultrafast surface plasmon (SP) dynamics in small noble-metal nanoparticles. 
Section ^ concludes the review. 
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2 Nonlinear Optical Dynamics in Modulation 
doped Quantum Wells 

By separating the ionized impurity donors from the electrons, modulation 
doping can produce a high mobility two-dimensional electron gas. In contrast 
to the case of a photoexcited electron gas, the temperature of such a system 
can be lowered, which enables the observation of fascinating many-body 



phenomena such as the Quantum Hall Effect ^ |33| and the Fermi 
Edge Singularity f^. 

Recently there has been a lot of interest in studying the ultrafast dy- 
namics in MDQW's. The role of the interactions between the photoexcited 
and the Fermi sea (FS) carriers during the thermalization stage II has been 
reviewed e.g. in Refs. P, fOl and is only briefly discussed here. In the 



section |3.2| we will review corresponding measurements in metal films and 
nanoparticles and compare with MDQW's. Much less is known about the 
role of many-body effects in the coherent regime (stage I). Even though some 
aspects of the effects of the e-e scattering have been described within the 
dephasing time approximation, the role of the e-h correlations between the 
photoexcited holes and the FS electrons has only recently been studied. The 
main difficulty in treating such e-h correlation effects comes from the non- 
perturbative nature of the Fermi edge singularity (FES), which dominates 
the absorption spectrum close to the absorption onset, and which cannot be 
described within the dephasing time approximation. Important here is that 
the e-h interactions between a heavy photoexcited hole and the FS-electrons 
lead to the scattering of a macroscopic number of low-energy FS-pairs, which 
readjusts the entire FS into a new orthogonal scattering state in the course 
of the optical excitation (Anderson orthogonality catastrophe Before 
we discuss the ultrafast dynamics due to such effects, we briefly summarize 
in the next section the linear optical properties of MDQW's. 

2.1 Fermi Edge Singularity in Linear Absorption 

Close to the onset of linear absorption and for low-temperatures, the optical 
properties of MDQW's are dominated by the FES. The latter is a many-body 
resonance that has been observed in doped semiconductors (see e.g p^, pO, 



38, 39, p|) as well as in metals [HOl Ell, where it is referred to as the X- 



ray Edge Singularity. Despite the screening of the Coulomb interaction, the 
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strength of this resonance is comparable to that of the undoped QW excitons. 
The non-Lorentzian hneshape of the FES can be viewed as originating from 
the decay of an excitonic bound state caused by its interactions with the 



a{ij)ocAf(—) , (4) 



gapless FS excitations In the case of holes localized due to e.g. the 

disorder or for core holes in metals, the lineshape of the FES close to 
the onset of absorption can be approximately described by using the analytic 
power law expression |^ 

where A/" is the density of states, Ep is the Fermi energy, uj is the optical 
frequency measured from the Fermi edge, and jS = 26/7T — is the 

FES exponent, where 6 ~ tan~^(7r5') is the s-wave phaseshift of the screened 
e-h potential V evaluated at Ep, and g = VAf is the dimensionless pa- 
rameter characterizing the e-h scattering strength. Such a non-Lorentzian 
lineshape results from the competition between the Mahan singularity, due 
to the attractive interaction between the FS and the localized hole (vertex 
correction effect, 0) and the Anderson orthogonality catastrophe due 
to the readjustment of the FS density profile during the optical transition 
(hole self-energy effect, 0]). In the case of finite hole mass rrih, the FES 
is broadened by an additional energy width of the order of the hole recoil 
energy ~ me/mhEp, where rrie is the electron mass ||4^, 

In a first approximation, the excitonic effects in a MDQW can be de- 
scribed by extending the mean field approach (HFA) to include the effects of 
the screening and the Pauli blocking due to the FS ^ . However, such 



a static treatment of the FS leads to a spurious bound state 0, ^ with re- 
spect to the Fermi energy, Ep, referred to in the following as the HFA bound 
state This discrete state, with binding energy Em, would appear at the 
energy Ep — Em- Obviously, for Em < Ep (as in typical MDQW's), such a 
state cannot exist since it overlaps with the FS continuum, with which it in- 
teracts via the e-h potential 0, ^ . The "unbinding" of this HFA bound state 
occurs via its interactions with the FS excitations, which are not taken into 



account in the HFA Note that, in two-dimensional systems, a static FS 
allows for bound states even for arbitrarily small attractive interactions and 
therefore such an unbinding cannot arise from static screening or from Pauli 
blocking effects. This spurious bound state could be artificially merged with 
the continuum by introducing a dephasing time comparable to its binding 
energy. Such an approximation, however, neglects competely the dynamical 
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correlations between the photoexcited e-h pair and the FS excitations. A 
microscopic description of the unbinding of the HFA bound state presents 
a nontrivial problem due to the non-perturbative nature of the e-h corre- 
lations between the photoexcited hole and the FS excitations |^2|]. Even in 
linear absorption, within Green function techniques, this problem is rather 
involved because vertex correction diagrams with arbitrarily many crossed 
e-h interaction lines are as divergent as the ladder diagrams and should be 
treated on equal footing [|^, To perform such a task, one must sum up at 
least the parquet diagrams and address the three-body correlations between 
the photoexcited hole and a FS excitation [^, 0. Therefore, alternative 
methods were developed for the case of linear absorption, based on Fermi's 
golden rule with many-electron eigenstates expressed in terms of Slater de- 
terminants [^. Such approaches become exact in the limit rrie/mh -C 1 and 



41, 50, 48, 51 



cluster expansion (CCE) |£2| 



describe quite accurately the FES lineshapes observed in typical MDQW's 
Another approach to the FES problem is based on the coupled 
|, 0, 55|. This general many-body technique 
provided exact results in the limit nie/mh <C 1 ||5^, and was used to treat 
the hole recoil correlations in one dimension ||5^. In the latter case, an exact 
solution was obtained for rrie = rrih- The CCE has also been used to describe 
the e 



-e correlation effects 60 



More importantly, however, this method 
is well-suited for describing correlations in non-equilibrium systems, where 
it retains the advantages of diagrammatic expansions without resorting to 



the Kadanoff-Baym ansatz or to the Markovian approximation pT 



2.2 Ultrafast Dynamics 

In this section we summarize some recent studies of the dephasing and re- 
laxation processes in MDQW's. The earlier work focussed on the role of the 
interactions between the FS and the photoexcited carriers on the thermal- 
ization during relaxation stage II. In most cases, a non-thermal distribution 
of real electron and hole carriers was photoexcited well above the Fermi 
surface and the subsequent incoherent population dynamics was monitored 
with ultrafast pump-probe spectroscopy. Note that, in such experiments, 
the measured differential transmission is 

n^TL . ^\ '^t.(u,,t) ue,) - ue, = 0) 
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where Ts{Ep) is the transmission coefficient in the probe direction in the 
presence of the pump field £p and r denotes the time delay between the 
pump and probe optical pulses. For samples with sufficiently small thickness 
d such that Aa d is small, the differential transmission reproduces the change, 
Aa(a;, r), in the probe absorption coefficient r) which is induced by the 
pump photoexcitation: DST{u!,t) ~ —Aa{uj,r)d. 

For a simple intuitive interpretation of the thermalization experiments in 
MDQW's, it is often assumed that, for free carriers photoexcited well above 
the onset of absorption and for time scales much longer than the dephas- 
ing time, the differential absorption spectrum maps the carrier distribution 
functions at the probe photon energy at time r. In particular, the following 
approximation is often used: 

Aa(cu,r) = (l-/e-A)aM (6) 

where a{uj) is the linear absorption coefficient and /e and fh are the electron 
and hole distribution functions at the corresponding energies. Within such 
an approximation, the time evolution of the differential absorption is deter- 
mined by the carrier distribution functions at time r. There have been many 
theoretical calculations of the distribution function time evolution, mostly 
based on numerical solutions of the semiclassical Boltzmann equations (see 



e.g. Refs. |10, 14, p|). A detailed review of such an approach to thermaliza- 



tion can be found in Ref. []TU[. More recently, the dynamics in the quantum 



kinetic regime and the shortcomings of the Boltzmann equations have been 
addressed, as reviewed e.g. in Ref. WU . 



Knox et. al. studied the effect of the inelastic e-e scattering on 

the thermalization by photoexciting a non-thermal carrier distribution well 
above the Fermi surface and then monitoring its time evolution to a thermal 
distribution at the bottom of the band. They observed that the e-e inter- 
actions between the photoexcited and FS electrons significantly enhance the 
thermalization rates as compared to undoped semiconductors. In particular, 
the thermalization time in MDQW's was found to be of the order of ~10fs. 
The differential transmission lineshape mimicked that of a Boltzmann dis- 
tribution peaked at the bottom of the conduction band, which indicates a 
non-degenerate carrier distribution. This is expected at room temperatures 
(as in the experiment of Ref. |^), where the thermal energy exceeds the 
Fermi energy, and at high intensities, where the photoexcited carrier den- 
sity exceeds the FS density. Knox et. al. also observed an instantaneous 
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redshift due to the bandgap renormalization by the photoexcited carriers 
0. However, as we discuss later, the situation changes drastically at low 
temperatures, where the FS electrons have a sharp Fermi-Dirac distribution. 
Wang et al. |Q performed low temperature pump-probe measurements and 
found that, at low excitation densities, the differential transmission lineshape 
corresponded to a redshift of the FES due to the incoherent bandgap renor- 
malization. Importantly, they showed that such a redshift slowly builds up 
with time as the hot electron temperature decreases due to the cooling to 
the lattice temperature. This observation points out that the bandgap renor- 
malization depends not only on the carrier density but also on the carrier 
distribution. At high photoexcited densities, the heating of the electron gas 
becomes stronger and Wang et. al. observed a bleaching of the FES due to 
the smearing of the sharp Fermi surface. 

Both the above experiments were interpreted in terms of e-e interactions 
between the photoexcited and the FS carriers. Tomita et. al. performed 



ultrafast luminescence measurements at low temperatures to investigate the 
role of the e-h interactions on the thermalization in n-type MDQW's. In 
particular, they measured the time evolution of the luminescence intensity at 
the bandgap energy, where the optically-induced perturbation of the electron 
distribution due to, e.g., heating is minimal. Tomita et. al. found that the 
number of holes that thermalized to the top of the valence band as a function 
of time could not accounted for if one assumes a thermalized hole distribution. 
This result indicates that a substantial fraction of the photoexcited holes are 
non-thermal for time intervals as long as ~ 800fs. Furthermore, Monte- 
Carlo simulations indicated that the e-h scattering is the dominant hole 
thermalization mechanism. Similar measurements in n-doped bulk GaAs 
were performed by Chebira et. al. (at low temperatures) and Zhou et. 
al. (at room temperature) [Q. Woerner et. al. ^ studied the hole 
relaxation in p-doped semiconductors by photoexciting holes from the heavy- 
hole to the split-off valence band and then monitoring the time evolution of 
the inter-valence-band absorption spectrum by using femtosecond pulses in 
the mid-infrared region. Due to the absence of conduction electrons, the 
dynamics is then solely determined by the h-h and h-p scattering. They 
found that the inter-valence-band scattering via the emission of e.g. optical 
phonons occurs within a short time interval <100fs, and that subsequently 
the phototexcited holes slowly thermalize via h-h scattering with the hole 
FS within a time interval ~ 700fs. 

An important effect of the quantum confinement in QW's is the formation 
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of discrete conduction and valence subbands. Recently, there has been much 
interest in exploring the role of the intersubband excitations in the ultra- 
fast dynamics, partly motivated by device applications such as the quantum 
cascade laser Lutgen et. al. [^, ^] performed femtosecond pump- 



probe measurements and showed that the intersubband nonlinear absorption 
is determined by both the intersubband and the intraband electron relax- 
ation. They measured the effects of intersubband pump excitation on the 
nonlinear absorption due to probe-induced interband transitions from the 
valence to the conduction band. They observed a population relaxation that 
proceeds in two stages. First, the electrons photoexcited by the pump from 
the Fermi sea in the lowest conduction subband relax back via intersubband 
scattering within a time interval ~ 1.3 ps. The latter is consistent with the 
intersubband relaxation times due to longitudinal optical phonon emission 
|]72| , [75[] . This process is followed by intraband e-e relaxation similar to the 
experiments discussed above. Lutgen et. al. also measured the time evolu- 
tion of the nonlinear absorption spectrum due to optical transitions between 
the conduction subbands. They observed an intial decay of the differential 
transmission within a time interval ~ 2-3 ps, which they attributed to the 
population relaxation due to the intersubband scattering. This was followed 
by a slowly decaying differential transmission signal, whose time evolution 
depended strongly on the photoexcitation frequency (within the linear ab- 
sorption linewidth). The latter regime was attributed to the intraband popu- 
lation relaxation within the lowest subband as the carriers cool to the lattice 
temperature. 

We now come to the effect of a cold FS on the e-e scattering time. One 
would expect that the latter should become short as the FS density increases. 
This is indeed the case for low FS densities, high temperatures, or for large 
excess energies from the Fermi surface. Hawrylak et. al. obtained 
the energy-dependent e-e scattering times by calculating the equilibrium 
self-energies including the effects of the electron-plasmon and e-p interac- 
tions. They predicted step-like decreases in the scattering time for electron 
energy above the Fermi surface that exceeds the onset of plasmon and opti- 
cal phonon emission. They also concluded that short-range e-e correlations 
(not included in the RPA) can significantly affect the scattering times at high 
energies. Such fast e-e scattering of photoexcited electrons well above the 
Fermi surface can explain for example the absence of spectral hole burning in 
the pump-probe spectra of MDQW's |]62|, ^J. However, the situation changes 



drastically for electron energies close to the Fermi surface and temperatures 
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smaller than the Fermi energy. Under such conditions, the e-e scattering is 
in fact suppressed in the presence of a degenerate FS. Indeed, a photoexcited 
electron lowers its energy by interacting with a FS electron while the latter 
scatters above the Fermi surface. For electrons close to the Fermi surface, 
the drastic reduction of the phase space available for scattering due to the 
Pauli blocking by the FS electrons as well as the screening leads to the strong 
suppression of the e-e scattering This is due to the sharp Fermi-Dirac 
distribution of the FS electrons at low temperatures. In fact, at zero tem- 
peratures, the e-e scattering time Te{E) becomes infinite right at the Fermi 
surface, E = Bp |1|, . According to Fermi liquid theory , 

where fc^T is the thermal energy. Note that the above result applies both 
to metals and MDQW's due to the similar values of r^. As the electron 
temperature increases, the smearing of the Fermi-Dirac distribution allows 
for the scattering of a photoexcited electron to states below Ep and therefore 
Te{E) decreases as compared to the zero temperature limit. Even though the 
e-e scattering rate initially increases with electron density for low densities, 
it decreases at higher densities when the FS temperature is low and the 
Fermi edge is sharp |75[. Therefore, one should expect delayed dephasing and 
thermalization processes as the electrons approach the Fermi surface (note 
that the above phase space restriction does no apply to the e-h scattering of 
the photoexcited hole with the FS electrons). 

The effect of a cold FS on the dephasing time of the interband optical 
polarization was investigated by Kim et al. |76| using transient FWM spec- 
troscopy 0. In such experiments, two optical pulses separated by a time 
delay r propagate along two different directions, ki and k2, and interfere 
within the sample, thus generating a nonlinear polarization. In FWM spec- 
troscopy, one measures the signal that emerges along the direction 2k2 — ki . 
In the case of optically-thin samples, the time-resolved FWM signal corre- 
sponds to the square of the amplitude of the nonlinear polarization along the 
above direction 0. The time-integrated FWM signal (TI-FWM) is given 
by the integral over all times of the time-resolved FWM signal and is a func- 
tion of the time delay r. To a first approximation, the FWM signal can be 
interpreted similar to atomic and molecular systems based on a two-level 
system Within such a model, in the case of homogeneous broaden- 

ing, the time-resolved FWM signal is emitted immediately after the second 
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pulse arrives and decays with a decay time T2/2, where T2 is the dephasing 
time. The TI-FWM signal also decays with a decay time T2/2 and van- 
ishes for negative time delays where pulse 1 comes after pulse 2. In the case 
of inhomogeneous broadening, the time-resolved FWM signal is delayed by 
a time interval r after the second pulse arrives and corresponds to a pho- 
ton echo |jl8|, ^ |1^. The TI-FWM signal decays with a decay time T2/4, 
which is a factor of 2 smaller than in the case of homogeneous broadening, 
and also vanishes for negative time delays. In undoped semiconductors, sig- 
nificant deviations from the above simple picture have been observed and 
have been attributed to Hartree-Fock exciton-exciton interactions and to 



exciton-exciton correlations |jT2[. For example, such effects lead to a strong 
TI-FWM signal for negative time delays as well as a delayed time-resolved 
FWM signal. For a detailed review of the exciton-exciton interaction effects 
in undoped semiconductors, see e.g. Ref. |T^. 

In MDQW's, Kim et. al. measured long dephasing times of a few pi- 
coseconds for energies within the frequency range of the FES (i.e. close to 
Ep). These decreased to less than lOOfs as the excitation frequency exceeded 
the FES peak by an energy comparable to the Fermi energy. The above 
result is consistent with the Fermi liquid energy dependence of the e-e scat- 
tering time, Eq. |^, and points out that, unlike for the e-h correlations, in 
MDQW's the e-e scattering is suppressed within the frequency range of the 
FES. This experiment also suggests that the hole dephasing times are of the 
order of a few picoseconds. In the time-resolved FWM signal, Kim et. al. 
observed photon-echo-like behavior, which indicates that, despite the strong 
FES peak in the linear absorption, the continuum of interband e-h states 
leads to behavior similar to that of an inhomogeneously-broadened system. 
Finally, in MDQW's, the pair-pair interactions, analogous to the exciton- 
exciton interactions in undoped semiconductors, are screened out and the 
negative time delay TI-FWM signal characteristic of exciton-exciton inter- 



actions was absent for high Fermi sea densities |]76|. The above experimental 



results were interpreted by Hawrylak et. al. |7j] within the two-level system 
approximation (which neglects all FES correlation effects) by introducing 
energy-dependent dephasing times determined by the e-e and e-p scattering 
times. 

In order to provide further insight into the dephasing processes in MDQW's, 
Bar-Ad et. al. performed low temperature FWM measurements in the 



presence of a magnetic field parallel to the QW growth axis. For zero mag- 
netic field, they observed a dephasing time ~ Ips. For magnetic fields such 
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that only the lowest Landau level is occupied, they observed a suppression of 
the e-e scattering. In particular, the dephasing times were then found to be 
longer, ~1.5 ps, and, unlike in the zero field case of Ref. [0, they did not 
depend significantly on the frequency of the optical excitation. Bar Ad et. 
al. attributed the above observations to the suppression of the phase space 
available for scattering induced by the magnetic confinement. Indeed, the 
application of a magnetic field changes the continuous density of states into 
a series of discrete highly degenerate Landau level peaks, broadened by the 
disorder. Finally, Bar-Ad et. al. observed a FES in the FWM spectrum at 
high magnetic fields, which they attributed to the optical transitions to the 
occupied lowest Landau level. 

The dephasing of the intersubband optical polarization in MDQW's was 
recently investigated by Kaindl et. al. using transient FWM in the 

midinfrared. They observed an exponential decay of the TI-FWM signal that 
was fairly insensitive to the temperature. By modelling their data with a HFA 
calculation, Kaindl at. al. deduced the values of the dephasing times to be of 
the order of hundreds of femtoseconds and concluded that the main dephasing 
mechanism is the e-e scattering. By increasing the Fermi sea density, they 
observed faster dephasing. In another experiment, Bonvalet et. al. 
used an extremely short, ~ 12fs, pulse to excite a wavepacket consisting of 
states in the two lowest conduction subbands of a QW. By measuring the 
radiated coherent electromagnetic field, they observed a quantum beat at 
room temperature oscillating at the inverse intersubband excitation energy. 
From the decay of the signal they deduced an intersubband polarization 
dephasing time of llOfs in the case of a MDQW, which was shorter than 
the dephasing time of 180fs observed in an undoped QW. Such polarization 
dephasing is determined by the destructive interference effects as well as by 
the different scattering processes. 

The above measurements of the polarization dephasing were performed 
under resonant photoexcitation conditions. In contrast to undoped semicon- 
ductors, the relaxation of such real pump-induced e-h pairs due to their e-e 
interactions with the FS will obscure any coherent effects. Fortunately, the 
dissipative processes can be suppressed by tuning the pump frequency below 
the onset of absorption. In this case, only virtual carriers are excited by 
the pump field and thus the coherent effects will dominate. For well-below- 
resonance pump excitation, the e-h pair phase is primarily determined by the 
optical field and the time dependence of the interband polarization follows 
that of the pump pulse (adiabatic following |19]). In this case, the pump 



18 



creates a truly coherent e-h pair many-body state that lasts for the duration 
of the pulse. Such pump-probe experiments in atomic systems showed a res- 
onance blueshift ( optical Stark effect) accompanied by a resonance bleaching 
due to the Phase Space Filling by the pump-induced e-h pairs [|19], |80 



In the case of excitons, it was shown that the Coulomb interaction can sig- 
nificantly alter the bleaching observed in atomic systems. In particular, at 
low pump intensities or for pulse duration longer than the dephasing time, 
the interactions lead to a pure exciton blueshift without significant bleach- 
ing. For high intensities and pulse duration shorter than the dephasing time, 
the exciton bleaching is very strong. The exciton ac-Stark effect has been 
reviewed e.g. in Refs. 0, |^, Such interaction effects in undoped 



semiconductors suggest that the many-body e-h correlations between the 
photoexcited holes and the FS electrons will significantly affect the coher- 
ent nonlinear optical response of the FES. Brener et. al. |S^, |S3I] studied 



this issue by performing pump-probe measurements for pump detunings ~ 
50meV below the Fermi surface, much larger than the e-h Coulomb energy 
Em (which is much smaller than the Fermi energy ~15meV). They observed 
a qualitatively different bleaching of the FES as compared to excitons in 
undoped QW's, which they attributed to the distinct nature of the two reso- 
nances. Schafer studied this issue by performing quantum kinetic calcu- 
lations using Keldysh Green's functions. Such an approach treats the HFA 
coherent effects and includes the scattering contributions within the second 
order Born, Markovian, and static screening approximations. The numeri- 
cal results for Fermi energy equal to the exciton binding energy showed a 
strikingly different bleaching between the FES and the exciton, which was 
attributed to the e-e scattering. 

As we discuss at length in section ^, the e-h correlations between the 
photoexcited hole and the FS electrons lead to strong dephasing and affect 
qualitatively the coherent ultrafast dynamics of the FES. An important dif- 
ference from undoped semiconductors is that, due to its gapless excitation 
spectrum, a FS responds unadiabatically to time-dependent perturbations. 
In contrast, because of its finite Coulomb binding energy, an exciton can be 
polarized by the pump optical field without being ionized. A non-equilibrium 
treatment beyond the HFA is necessary in order to take into account the una- 
diabatic time-dependent change in the e-h pair-FS interactions and the e-h 
scattering processes induced by the ultrafast pump excitation in the coherent 
regime. As has already been noted in the case of photoexcitation of an elec- 



tron gas within the continuum of states of an undoped semiconductor |84|, 
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the loss of coherence due to the many-body correlations cannot be fully de- 
scribed within the dephasing time approximation. In section ^ we show that 
memory effects due to the e-h correlations lead to a time evolution of the 
pump-probe signal characterized by the inverse Coulomb energy Em rather 
than by the dephasing time (as in the case of exciton bound states). 

Concluding this section, let us briefly discuss a very recent FWM exper- 
iment by Fromer et. al. which demonstrated for the first time non- 
Markovian memory effects in the Quantum Hall Effect regime. Even though 
for magnetic fields between 5.5 and 6.5 T the TI-FWM profile was found to 
be a single exponential with an unusually long decay time, for magnetic fields 
that exceed 7.5 T such a profile was more complicated and characterized by a 
change of slope indicating memory effects in the polarization dynamics. Such 
effects were also seen in the frequency domain: the FWM spectrum profile 
changed from a Lorentzian lineshape to an asymmetric one corresponding to 
a frequency dependent width, r(cc;). To interpret such behavior, we note that, 
at high magnetic fields such that the cyclotron energy, huc-, is large compared 
to other characteristic energies of the system, the relaxation is dominated by 
intra-Landau-level processes. Such scattering by collective excitations in- 
volves the matrix elements of the dynamically screened interaction, U^{t, t'), 
which in the lowest Landau level have the form: 

U<{t,t') = /(0^e-^'^/\^x<(t,Oc.,(g), (8) 

where Xq{t,t') = (pq(t')p_q(t)) is the density-density correlation function 
projected onto the lowest Landau level [^, and Pq(t) is the corre- 
sponding density operator. Here, Vg is the unscreened Coulomb interac- 
tion, / = [fi/eBy/'^ is the magnetic length, and the coefficients Cjj(g) with 
z,j —>■ e,h model the asymmetry in the e-e and e-h interaction matrix el- 
ements, which originates from the difference between the electron and hole 
lowest Landau level wavefunctions. Because of the breakdown of perturba- 
tion theory due to the Landau level degeneracy in 2D systems, it is incorrect 
to evaluate x^{tyt') within the standard RPA []TU[. Instead, one should ac- 



count for the true excitations of the interacting two-dimensional electron 
liquid. Several models can be found in the literature, and we base our dis- 
cussion on the magnetoroton model, which is the one best suited for the 
filling factors u in the experiment of Ref. [^. The most salient features 
are, however, general and model independent. The magnetoroton dephasing 
mechanism is somewhat similar to that of acoustic phonon scattering. Under 
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the experimental conditions of Ref . , to a very good approximation, the 



intra-Landau-level collective excitations are not affected by the small density 
of photogenerated carriers, so one can use the equilibrium density correlation 
function[^. The equations for the density matrix elements then read 



dpij 



= I 

scatt k 



y: f dt'Gi{t-t')G'^it'-t) 

1 J —oo 



dt 

X {[U<,{t - t') - U<{t - t')]p<{t')p>^{t') - (<->) j , (9) 

where G\^"'{t) is the retarded/advanced Green function, pf^ = pij, and pfj = 
6ij — pij. If all Uij are equal, i.e., Cjj(g) = 1, then the polarization scattering 
term vanishes This corresponds to identical electron and hole wavefunc- 
tions in the lowest Landau level. In practice, there is always an asymmetry 
between electrons and holes, due to, e.g., differing band offsets, lateral con- 
finement, and disorder. Using the results of Ref. |p6|) Eq. (pi) takes the form 



U<(t) = f ^S_e-9''V2^2 . ) 



xsq[{Nq + l)e'^^' + Nqe-'^^% (10) 

where Nq is the Bose distribution function for magnetorotons of energy ujq, 
and Sq is the static stucture factor of the 2D electron liquid in the lowest 
Landau level. By comparing Eqs. (|]) and (|TOl), we see that the u dependence 
of r(cc;) is determined by the Fourier transform of U^{t), which is governed 
by the q dependence of Sq. In the lowest Landau level, we have Sq = (2z/~^ — 
l)sq, where Sq ~ (ql)"^ for g/ -C 1, ~ exp(— g^/^/2) for ql ^ 1, and Sq 
displays a peak for g/ ~ 1 that leads to the magnetoroton excitations. The 
corresponding resonance in r(ci;) near the magnetoroton energy leads to non- 
Markovian behavior with a characteristic response time of approximately 
the inverse of this energy. The latter is estimated from the gap A at the 
magnetoroton dispersion minimum, A ~ 0.1(e^/e/) for the range of u in this 
experiment which for i? = 10Tis^l.5 meV. 

This concludes our overview of the main ultrafast dynamical features 
observed in MDQW's. In the next section we discuss the ultrafast dynamics 
in the case of metal nanoparticles and compare to MDQW's. 
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3 Nonlinear Optical Dynamics in Metal Nanopar- 
ticles 



The properties of small metal particles in the intermediate regime between 
bulk-like and molecular behavior have been the subject of great interest 
lately. This was motivated in part by the need to understand how the prop- 
erties of matter evolve at the transition from atoms and molecules to bulk 
solids. An additional motivation comes from the technological trend towards 
electronic and optoelectronic devices based on smaller and smaller solid state 
structures |5D|, ^ Metal clusters are also being used in a variety 

of applications, ranging from catalysis to biological and medical applications 
( see e.g. Refs. |]^, 0, ^]). The electronic and thermodynamic 

properties of metal clusters have been extensively reviewed, e.g. in Refs. 
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lOOj , |101| , p.02|| . The linear optical properties of metal nanoparticles have 
also been reviewed in detail by Kreibig and Vollmer ||103| . 

It has been known for a long time that surface collective excitations play 
an important role in the absorption of light by metal nanoparticles. In large 
particles with sizes comparable to the wave-length of light A (but smaller 
than the bulk mean free path), the lineshape of the surface plasmon (SP) 
resonance is determined by the electromagnetic effects |p.04] , |105| , |106| , |103|| . 
However, as the size of the nanoparticle becomes smaller than the mean free 
path of electrons in the bulk metal, quantum confinement becomes impor- 
tant. In small nanoparticles with radii R <^ X, the absorption spectrum is 
governed by quantum confinement effects. For example, the momentum non- 
conservation due to the confining potential leads to the Landau damping of 
the SP and to a resonance linewidth inversely proportional to the nanoparti- 
cle size 
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A review of the extensive theoretical and experimental studies 
of this effect may be found in Ref . [ p.03|| . More recently, it was established 
that the static nonlinear optical properties of small nanoparticles are also 
affected by the confinement. In particular, a size-dependent enhancement of 
the third-order nonlinear optical susceptibilities for monochromatic photoex- 
citation, caused by the elastic surface scattering of single-particle excitations, 
was predicted by Flytzanis and collaborators ||108| , |109|| and observed exper- 
imentally by Yang et. al. ||110|| . Dielectric confinement also enhances the 
optical nonlinearities close to the SP frequency ||109| , |103| , Before we 

proceed with the discussion of the ultrafast dynamics of such confined Fermi 
seas, let us briefly summarize in the next section the main features of the SP 
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resonance in the linear absorption spectrum. 



3.1 Surface Plasmon Resonance in Linear absorption 

In this section we summarize the basic facts regarding the hnear absorption 
by small metal particles embedded in a medium with dielectric constant e^- 
We will focus primarily on noble metal particles containing several hundreds 
of atoms; in this case, the confinement affects the extended electronic states 
even though the bulk lattice structure has been established. When the parti- 
cle radii are small, i? <^ A, so that only dipole surface modes can be optically 
excited and non-local effects can be neglected, the optical properties of this 
system are determined by the dielectric function ||103|| 

ecoi(u;) = + 3pe™ ^(^) ^ (^l) 
e{uj) + 2em 

where e{uj) = e'{uj) + ie"{uj) is the dielectric function of a metal particle 
and p <^ 1 is the volume fraction occupied by nanoparticles in the colloid. 
Since the d-electrons play an important role in the optical properties of noble 
metals, the dielectric function e{uj) includes also the interband contribution 
ed{uj) due to transitions from the d-band to the s-p conduction band. For 
p <C 1, the absorption coefficient of such a system is proportional to that of 
a single particle and is given by ||103 | 



a{Lu) = -9pe'J,^-lm^, (12) 

where 

€s{uj) =ediuj)-ujl/uj{uj + i'js) + 2em (13) 

plays the role of an effective dielectric function of a particle in the medium. 
Its zero, e'g{ujs) = 0, determines the frequency of the SP, Ug. In Eq. (|13|), 
ujp is the bulk plasmon frequency of the conduction electrons, and the width 
7s characterizes the SP damping. The semiclassical result Eqs. ( |12|) and 
([T3| ) apphes to nanoparticles with radii R ^ q~^, where g^,^ is the Thomas- 
Fermi screening wave-vector {q~^ ~ 1 A in noble metals). In this case, the 
electron density deviates from its classical shape only within a surface layer 
occupying a small fraction of the total volume ||112|| . Quantum mechani- 



cal corrections, arising from the discrete energy spectrum, lead to a width 



7s ~ "^f/-^' where = kp/m is the Fermi velocity [103, 107]. Even though 
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Isl^s ~ {<1tf^)~^ ^ 1; this damping mechanism dominates for sizes i? < 10 
nm. Additional contributions to such a width come from the e-p, e-e, and 
electron-impurity interactions as well as from the disorder. On the other 
hand, in small clusters containing several dozens of atoms, this semiclassi- 
cal approximation breaks down and density functional or ab initio methods 



should be used 103, 99, 100, 101, |102|, 113|. In particular, discrete electronic 



levels whose width does not exceed their spacing are expected for particle 
sizes smaller than a few angstroms. In the latter regime, the e-e interactions 
are similar to those in atoms and molecules. One therefore expects that the 
role of the many-body correlations increases as we approach the crossover 
from quasi-continuous to discrete nanoparticle energy levels. 

It should be noted that, in contrast to the scattering with surface collec- 
tive excitations, the e-e scattering is not too sensitive to the nanoparticle 
size as long as the condition g^^i? ^ 1 holds ||114|| . Indeed, for such sizes. 



the static screening is essentially bulk-like. At the same time, the energy de- 
pendence of the bulk e-e scattering rate (Eq. ^ 7e oc {E — EpY comes from 
the phase-space restriction due to the energy and momentum conservation, 
and involves the exchange of typical momenta q ~ q^^. If the size-induced 
momentum uncertainty 6q ~ R^^ is much smaller than q^^, the e-e scat- 
tering rate in a nanoparticle is not significantly affected by the confinement 



11511 . Below we will see however that this is not the case for the electron-SP 



interactions. 



3.2 Ultrafast Dynamics 

Even though the electronic, thermodynamic, and optical properties of metal 
nanoparticles have been extensively studied, the role of confinement in the 
electron dynamics is much less understood. Examples of outstanding issues 
include the role of e-e interactions in the process of cluster fragmentation, 
the role of surface lattice modes in providing additional channels for intra- 
molecular energy relaxation, the influence of the electron and nuclear motion 
on the superparamagnetic properties of clusters, and the effect of confine- 
ment on the nonlinear optical properties and transient response under ultra- 
fast excitation ||103| , |113| . These and other time-dependent phenomena can 



be studied with femtosecond nonlinear optical spectroscopy, which in these 
structures provides time resolution shorter than the relaxation times. Similar 
experimental studies in bulk metals and metal interfaces and surfaces have 



been reviewed e.g. in Ref. ^Tj. Here we will focus on the more recent work 
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in metal nanoparticles. 

Extensive experimental studies of the electron relaxation in noble-metal 
nanoparticles have recently been performed using ultrafast pump-probe spec- 
troscopy. In contrast to the situation in semiconductors, the dephasing pro- 
cesses in metals are very fast. In metal nanoparticles, the dephasing time 
can be deduced, e.g., from the linear absorption SP width to be of the order 
of a few fs. Such measurements were performed by Klar et. al. ||116|| , who 
distinguished between homogeneous and inhomogeneous broadening by mea- 
suring directly the SP lineshape of a single nanoparticle using a near-field 
antenna effect. The decay time of the SP resonance has also been studied 
with second and third harmonic generation measurements | 117 , 118| . 

Within several femtoseconds after the photoexcitation, the optically- 
excited electrons and holes form non-equilibrium populations. The pump- 
induced intraband transitions create a broad non-thermal electron distribu- 
tion that extends from the Fermi energy up to the pump photon energy, 
while the interband transitions create additional electron and d-band hole 
populations when the pump frequency exceeds the d-band to conduction 
band transition threshold. Within a few fs, the high energy electrons scatter 
to lower energies due to e-e interactions with the FS electrons. The latter 
interactions are screened within time intervals of the order of the inverse 
plasma frequency, typically ~ Ifs. Note that, similar to MDQW's [0, at 
high electron excess energies the e-e scattering time is of the order of a fs, 
much shorter than the e-p scattering times (of the order of a ps). However, 
as such electrons scatter down to the Fermi surface, the e-e scattering times 
become much longer, of the order of several hundreds of fs. As discussed in 
the previous section, this is a consequence of the Pauli blocking of the phase 
space available for scattering and the screening of the e-e interaction, which 
leads to energy-dependent scattering times given by Eq. |^ |jT|. Therefore, 
one expects that thermalization should slow down as the electrons approach 
the Fermi surface. In the case of metal films, a non-thermal population con- 
sisting of a hot Fermi-Dirac distribution together with a tail of high energy 
electrons right above the Fermi surface (and a tail of low energy holes be- 
low the Fermi surface) was indeed observed in time-resolved photoemission 
experiments during time intervals of the order of hundreds of fs ||119| , |120| . 
Pump-probe experiments also showed a clear signature of delayed thermal- 
ization (see e.g. ||121| , |122| , |123| , |124| , |125|| ). In particular, the rise time of 



the differential transmission and reflectance signal was observed to be of the 
order of hundreds of fs, much longer than the pulse duration. 
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A similar slow rise time of the differential transmission was also observed 
in the case of copper ||126|| , gold ||127|| , and silver [ p.28|| nanoparticles as well as 
in gold nanoshells ||129| , |13Q| . It has been reported that, in small nanoparticles 
with diameters ~5 nm, such a delayed thermalization regime (stage II) lasts 
longer than in metal films ||126| , |128|| . At the same time, the e-p interaction 
for such sizes was found to be weaker as compared to the bulk [|126|, |128|| . A 



possible explanation of the reduced e-p interaction is that, for the small sizes, 
the characteristic phonon energy, given by the Debye frequency, becomes 
smaller than the spacing between the nanoparticle energy levels close to the 
Fermi surface, which is of the order of vp/R- As a result, the scattering 
of an electron via bulk-like phonons is suppressed |131|. As the electrons 



approach the Fermi surface, the e-e scattering times become comparable to 
the e-p scattering times and therefore the e-p interaction also contributes to 
the thermalization of the electron gas. The suppression of the latter could 
lead to delayed internal thermalization of the electron gas. 

The main spectral feature in the differential absorption of the copper and 
gold nanoparticles (where the onset of interband transitions is very close to 
SP frequency) was a transient asymmetric broadening of the SP resonance 
g, |13|, |T3|, |T3|, |T3|, |126|, |127[ |133, g). Perner et. 



see e.g. 



Refs. 



al 



127|| observed a build-up of such a pump-induced broadening during 
time intervals of the order of Ips. Bigot et. al. ||126|| observed a similar 
buildup and also pointed out that the differential absorption lineshape could 
not be understood without including the effects of the energy-dependent e- 
e scattering on the interband dielectric function. They attributed the SP 
broadening mainly to the smearing of the Fermi-Dirac electron distribution 
close to the Fermi surface, due to the heating of the electron gas, which 
affects the interband dielectric function. Inouye et. al. ||137|| arrived at similar 



conclusions and argued that, in the case of resonant interband excitations, 
the contribution of the intraband e-e scattering to the damping of the SP, 
which is described by the width 7^ in the intraband dielectric function in Eq. 
([T3| ) ||125| , |139| , |140|| , plays a minor role. 

The time-delay of the SP broadening can be understood as follows. For 
quasi-equilibrium conditions, the nonlinear absorption spectrum can be de- 
scribed using the linear absorption results but with dielectric constants de- 
termined by the time-dependent carrier populations. As discussed above, 
within the time resolution of the experiment, the electron distribution de- 
viates from the equilibrium Fermi-Dirac distribution in the vicinity of the 
Fermi surface, due to the high electron temperature and the non-thermal 
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population [TT^, TW., 121 , 122|| . As can be seen from the linear absorption 
expressions, such a smearing of the FS distribution leads to the broaden- 
ing of the SP resonance. Therefore, the initial rise time of the differential 
transmission signal is determined by the non-equilibrium electrons ouside 
the equilibrium Fermi-Dirac distribution. Initially, such electrons occupy 
the broad photoexcited distribution and therefore their number is relatively 
small, determined by the pump intensity. However, their number increases 
with time as the e-e interaction leads to the scattering of more electrons 
out of the FS. This is the origin of the time-dependent increase in the SP 
broadening. 

Since the electron heat capacity is much smaller than that of the lattice, 
an electron temperature much higher than that of the lattice can be reached 
during subpicosecond time scales. Subsequently, the electron and phonon 
baths equilibrate through the e-p interactions over time intervals of a few 
picoseconds. During this incoherent stage, the hot electrons can be charac- 
terized by a thermalized Fermi-Dirac distribution with time-dependent tem- 
perature T{t), while the phonons can be characterized by a Bose-Einstein 
distribution with time-dependent temperature Ti{t). During the cooling of 



the hot electron gas, the SP width decreases ||132| , |133| , |134| , |135| , |136| , p.26 



1271 , |137| , |138|| . As soon as the electrons have equilibrated among themselves. 



one can study the subsequent time evolution by using a set of coupled dif- 
ferential equations for the temperatures T(t) and Ti{t), referred to as the 
two-temperature model |IiT| , P^ , |n7[: 



Ce(T)^ = V ■ («:eVT) - G(T - Ti) + P(r, t), 
at 

Q^=G{T-Ti), (14) 

where Ce{T) = 7T and Ci are the electron and lattice heat capacities, re- 
spectively, G is the energy transfer coefficient between the electrons and the 
lattice, which is proportional to the e-p coupling constant ||143|| , Kg is the 



electronic thermal conductivity, and the source term P(r, t) describes the 
local energy density per unit time absorbed by the electron system from the 
pump optical pulse. The first term on the rhs of Eq. (|T^) describes the ther- 
mal diffusion from the nanoparticle to the surrounding matrix, which occurs 
on a time scale of several tens of ps. The second term on the rhs of Eq. 



( [T^ ) determines the tranfer of heat from the electron gas to the lattice as the 



former cools down during a time interval of several ps. The above model as- 
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sumes electron and lattice temperatures larger than the Debye temperature 



142 1 and that the interactions can maintain a quasi-equilibrium for both the 



electron and the phonon populations at all times. The cooling of the electron 
gas to the lattice manifests itself via an exponential decay of the differential 
transmission. According to Eq. (|l^ , the decay time of the differential tran- 
mission signal is inversely proportional to the strength of the e-p interaction 



123| , |124|| . Stella et. al. [p.48|| performed transient reflectivity measurements 



in metallic tin nanoparticles with radii ranging from 2 nm to 6nm and saw 
a decay of the signal as a function of pump-probe delay that became faster 
by decreasing the nanoparticle radius. They deduced from their data a con- 
tribution to the energy relaxation time during the thermalization of the hot 
electron Fermi sea with the lattice that was inversely proportional to the 
nanoparticle radius. In a subsequent paper, Nisoli et. al. ||149|| performed 
femtosecond pump-probe measurements in solid and liquid gallium nanopar- 
ticles with radii ranging from 5 nm to 9 nm and observed energy relaxation 
time constants that varied from 1.6 ps to 600 fs with decreasing nanoparticle 
size. They observed similar electron relaxation dynamics in the solid and liq- 
uid nanoparticle phases, which indicates that in this system, the scattering of 
electrons with bulk phonons plays only a minor role in the relaxation because 
of the reduction in the available phase space due to the quantum confine- 
ment. Instead, the hot electrons transfer their excess energy to the lattice 
through the generation of surface vibrational waves (capillary waves), which 
leads to an electron-surface phonon interaction G inversely proportional to 
the nanoparticle radius ||150|| . An additional contribution to this effect may 



come from the fact that the heat diffusion to the matrix becomes faster with 
decreasing nanoparticle size due to a thermal buildup that occurs for larger 
particles. A similar effect was observed by Halte et. al. in the case of very 
small silver nanoparticles [T^]. The role of surface effects on the dynamics 
was also discussed in Ref. ||136| . 

In the gold and copper nanoparticles, condidered above, the SP frequency 
is very close to the onset of interband transitions. In contrast, in silver 
nanoparticles, the SP and interband transitions are well separated in energy. 
Halte et. al. ||128|| compared the electron dynamics between silver thin films 
and silver nanoparticles embedded in glass in the same spectral range, and 
studied its dependence on the photoexcitation intensity. In thin films, they 
found that the RPA dielectric function corresponding to the instantaneous 
temperature determined by the two-temperature model reproduced their re- 
sults quite well. In the case of the silver nanoparticles, the main feature 
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in the differential absorption spectrum was an apparent SP redshift, which 
was attributed to the changes in the real part of the interband dielectric 
function due to the thermal broadening of the electron distribution close to 
the Fermi level. It was also shown that the finite e-e lifetime of the carri- 
ers photoexcited via interband transitions was essential for obtaining such 
a redshift. Finally, the strong non-parabolicity of the conduction band at 
the energies corresponding to the optical frequencies was shown to be impor- 
tant for interpreting the data. Averitt et. al. ||129| , |130|| performed ultrafast 



pump-probe measurements in gold nanoshell structures, which consist of a 
dielectric core surrounded by a thin metallic shell of nanometer dimensions. 
In such structures, the SP frequency can be tuned by changing the ratio of 
the core diameter to the shell thickness, and the experiment was performed at 
frequencies far off the onset of the interband transitions. The delayed pump- 
induced broadening of the SP resonance suggests that a nonthermal electron 
population is present in these structures during the initial ~100 fs following 
the pump photoexcitation. The decay of the pump-probe signal was found to 
be somewhat slower than in the bulk, which was attributed to a weaker e-p 
interaction due to the reduction of phase space available for scattering with 
bulk phonons induced by the confinement. Finally, they observed a tran- 
sient blueshift of the SP resonance, attributed to the off-resonant interband 
transitions. 

Recent experimental results in the case of small noble metal particles 
126|, |128|| indicate that many-body correlation effects play an important role 



during the cooling of the electron gas to the lattice (relaxation stage III). 
Despite the similarities to the bulk-like behavior, observed, e.g., in metal 
films, certain aspects of the optical dynamics in nanoparticles are significantly 
different ||126| , |152| , |153|| . For example, the experimental studies of copper 
nanoparticles by Bigot et. al. [ p,26|| revealed that, for sizes smaller than ~ 
5nm, the decay times of the pump-probe signal depend strongly on the probe 
frequency in the immediate vicinity of the SP resonance. In particular, the 
relaxation is considerably slower at the SP resonance, and becomes faster 
right above and right below the SP frequency | ]126| , |152| , |153| |. This and other 
observations suggest that collective surface excitations play an important role 
in the electron dynamics in small metal particles. This important issue will 
be discussed at length in section 0. 
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4 Coherent Ultrafast Response of the Fermi 
Edge Singularity: Formahsm 

4.1 Basic equations 

Before we proceed with our discussion of the time evolution of the FES in 
the coherent regime, we outhne the main points of a formahsm |p.54| , p.55| , 
^ 156 , 15711 recently developed to account for the non-Markovian dynamics 



due to the dynamical FS response {e-h correlations). More details may be 
found in Appendix A. 

As discussed in section p.l| , even in linear absorption, the HFA and the 
dephasing time approximation have serious shortcomings when used to de- 
scribe the FES. In the case of a photoexcited electron gas, the correlation 
effects have been treated within the second Born approximation by using 
Keldysh Green functions (see e.g. Ref. [0]). However, at low temperatures 
and within the frequency range of the FES, the effects of the e-h interac- 
tion must be consistently accounted for to arbitrary order. For example, as 
discussed in section p.l| , even in linear absorption the vertex correction di- 
agrams with arbitrarily many crossed e-h interaction lines are as divergent 
as the HFA ladder diagrams and should therefore be treated on equal foot- 
ing. With Green functions, this requires summing up at least the parquet 



diagrams |^2|, a formidable task especially in the non-equilibrium fem- 
tosecond regime. Alternative methods are therefore desirable. In undoped 
semiconductors, the "dynamic controlled truncation scheme" has been used 
to treat the exciton-exciton correlations [11, ^ (see section 0). However, 



such a hierarchy of density matrix equations no longer truncates if the ground 
state of the semiconductor includes a FS. Furthermore, as we demonstrate 



in section |5.1j , an expansion in terms of the optical fields breaks down for fre- 
quencies within the FES range. Finally, in view of the significant complexity 
of the problem, it is highly desirable to use a method that also allows for a 
physically intuitive interpretation of the results. The purpose of this section 
is to outline the main points of such a method, which will be used in section 
^to describe the coherent nonlinear response of the FES. 

We first consider the case of one conduction subband and later extend 
the formalism to include the Coulomb coupling to a second subband. We 



also consider spinless electrons for simplicity. In the rotating frame [19|, the 
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Hamiltonian describing this system is 



H,,,{t) = H + Hp{t) + H,{t). 



(15) 



The first term is the "bare" Hamiltonian, 



H = Y. <^-q^-q + E(^q + ^Wcf^c, + ^ee + Vhh + Veh, 



(16) 



where is the creation operator of a conduction electron with energy 
and mass me, feLq is the creation operator of a valence hole with energy 
and mass m/j, Vee, V^/i, and Vhh describe the e-e, e-/i, and interactions, 



respectively, and Vt = Eg + Ef{1 + nie/mh) 



is the detuning of the 



central frequency of the optical fields, ujp, from the Fermi level. Eg being the 
bandgap (we set h = 1 everywhere). The second and third terms describe 
the coupling of the pump optical field, £p{t)e^^'''^~^'^''*, and the probe optical 
field, Ss{t)e'^''''~^'^p^^~'^\ respectively: 



Hp{t) = -i2Sp{t) [e'^^-^'U^ + H.c 
Hs{t) = -fiSsit) [e^'^-'-+*^*'"[/^ + H.c. 



(17) 



where the pump amplitude Sp{t) is centered at time t = and the probe 
amplitude £s{t) is centered at the time delay t = t, n is the interband 
transition matrix element, and 



^^=E*'-q 



(18) 



is the optical transition operator. The conventions for the time delay r are 
clarified in Appendix B. 

In many experiments, the amplitude of the probe field is much smaller 
than that of the pump, |i^p(t)| ^ |^s('^)|- In that case, as was shown in 
Ref. [154] (see Appendix A), the experimentally-measurable nonlinear op- 



tical polarization can be obtained in terms of the linear response functions 
of a "pump-dressed" semiconductor to a probe field (note that, within x^^\ 
this is true even for comparable pulse amplitudes). This "dressed" system is 
described by a time-dependent effective Hamiltonian H{t), which is obtained 
by performing a time-dependent Schrieffer- Wolff /Van Vleck canonical trans- 
formation on the Hamiltonian H + Hp{t) [[158| , |159| , |160| , |161|| . As we shall 
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see later, in all the cases of interest, the effective Hamiltonian H{t) has the 
same operator form as the bare Hamiltonian H, with the important differ- 
ence that the band dispersions (effective masses) and interaction potentials 
are time-dependent through Sp{t). Thus, the calculation of the nonlinear 
absorption spectrum reduces to that of the linear absorption spectrum of the 
"pump-dressed" semiconductor with uncoupled "effective bands" — a great 
simplification that allows us to use straightforward generalizations of well es- 
tablished theoretical tools in order to treat the correlations. In fact, such an 
approach mimics nicely the spirit of the pump-probe experiments and allows 
for a physically intuitive interpetation of the results. It is important to note 
here that this "pump-dressed semiconductor" approach is not restricted to 
monochromatic pulses and is valid for any pulse duration [|154 . 



The pump-probe nonlinear polarization has the following form (see Ap- 
pendix A): 



-iijp (t—r) 



dt'Ss{t'){Mt)\m}C{tX)u\t')\<i>oit')) ■ (19) 



Here, |<l'o(^)) is the state evolved with H{t) from the semiconductor ground 
state |0) of H, W{t) is the effective optical transition operator describing the 
probability amplitude for the photoexcitation of an e-h pair by the probe 
field in the presence of the pump excitation, and ]C{t, t') is the time-evolution 
operator satisfying the Schrodinger equation 



z^/C(t,t')=^W/C(t,0 



(20) 



The above equation describes the time evolution of a probe-photoexcited e-h 
pair in the presence of the pump excitation. The effective Hamiltonian and 
effective transition operator are given by (see Appendix A) 



H{t) =Ho + ^ {Sp{t) [V{t), U^] + H.c.) 



and 



+ 



1 r 



V\t),\u\V{t) 



(21) 



(22) 



where the operator V^{t), which generates the canonical transformation, sat- 
isfies the equation 



.dVHt)_ 
' dt 



H,V\t)] +ii£p{t)U\ 



(23) 
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with the initial condition = before the pump arrives. Eqs. pi| ) and 
([2^ ) include all the pump-induced corrections to f/^(t) and H{t) up to the 
second order in the pump optical field and are valid when (fiEp/Q) < 1 (for 
off-resonant excitation) or [fiEptp)"^ < 1 (for resonant excitation), where tp 
is the pump duration. 

It should be emphasized that, although Eq. ( pT]) gives the effective Hamil- 
tonian up to the second order in the pump field, SpS*, the polarization ex- 
pression Eq. (p!9| ) describes the effects of H in all orders. For example, as we 
shall see in section [4.3| , the pump-induced term in H contains self-energy cor- 
rections to the electron/hole energies, which describe (among other effects) 
the resonance blueshift due to the ac-Stark effect. Although the magnitude 
of these self-energy corrections, calculated using Eq. ([2T|), is quadratic in £p, 
the correct position of the resonance can only be obtained by evaluating the 
PP polarization (|T^) nonperturbatively (beyond x*'^'*)? i-G-? without resorting 
to the expansion of the time-evolution operator }C{t,t') in the pump field. 
Importantly, as we demonstrate in |5.1| , the same is true when calculating the 
effects of the self-energy corrections on the e-h correlations. As we shall see 
in section ^ such a nonperturbative (in the pump field) treatment of the non- 
linear response of the FES is crucial for the adequate description of the PP 
spectrum at negative time delays. In Sections 4^ and iA, we will describe 
the corresponding procedure, which accounts for the FS dynamical response. 
In contrast, the third-order polarization, x^^\ can be simply obtained from 
Eq. (|T^ by expanding )C{t,t') to the first order in the pump-induced term 
in H [second term in Eq. (pT])]. We did not include in Eq. (|1^) the biexci- 
tonic contribution (coming from the excitation of two e-h pairs by the pump 
and the probe pulses) since it vanishes for the negative time delays (r < 0) 
considered below where the coherent effects dominate [[154| • 

In addition to including important contributions beyond x^^^ via the solu- 
tion of Eq. (|20|) as discussed above, the advantage of Eq. (|l^), as compared 
to the equations of motion for the polarization, comes from its similarity to 
the linear polarization that determines the linear absorption spectrum |^ 



This can be seen by setting Sp{t) = in the Eqs. (^Tj) and (P^, in which case 
the effective time-evolution and optical transition operators transform into 
their "bare" counterparts: /C(t,t') — *• e"*-^*^*"* ^ and U^{t) W . Moreover, 
like W , the effective transition operator W{t) creates a single e-/ipair, while, 
as we shall see in section |4.3| , the effective Hamiltonian H{t) can be cast in 
a form similar to H. This allows one to interpret the Fourier transform of 
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Eq. (|T^ as the linear absorption spectrum of a "pump-dressed" semiconduc- 
tor with two uncoupled but time-dependent effective bands. This mapping 
simplifies significantly the calculation of the FES ultrafast nonlinear optical 
response by allowing a straightforward generalization of the CCE. It also 
allows one to interpret the various dynamical features in the nonlinear ab- 
sorption spectra, originating from the correlation effects, within the familiar 
framework developed for linear spectroscopy. 



4.2 Overview of the Coupled Cluster Expansion 

In this section, we show how the time-dependent CCE can be used to study 
the effects of the e-h correlations (dynamical FS response) on the time evo- 
lution of the e-h pair photoexcited by the probe. Our goal is to evaluate the 
many-body state = IC{t,t')W{t')\^Q{t')) that enters in Eq. This 

state satisfies the Schrodinger equation 

As already mentioned, H{t) has the same form as the bare Hamiltonian 
H. This allows us to obtain {"^{t)) through a straightforward generalization 



of the linear absorption calculations [55, |56|, ^ 58 1. After eliminating the 



valence hole degrees of freedom ||162| , |58| , is expressed in the CCE form 
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=e^W|$(t)), (25) 
where the time-dependent operator S{t) creates FS-pairs and is given by 

S{t) = s(p, k,t)aj,ak 

p>kp,k<kp 

+ ■S2(pi,P2,ki,k2,t)aJ,^aJ,^ak2aki H , (26) 

pi,P2>kp ,ki,k2<kp 



while the state |$(t)), discussed in Section describes the time evolution 
of the probe-induced e-h pair. In Eq. (p6D , the amplitude s(p, k, t) describes 
the e-h correlations which, in particular, are responsible for the unbinding of 
the HFA bound state; the two-pair amplitude S2 describes the e-e interaction 
effects at the RPA level and beyond. /^From a physical point of view, the 
operator c^^*^ describes the readjustment of the FS density profile during the 
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optical transition in response to the FS interactions with the photoexcited 
e-h pair. 

Substituting Eq. ( |25D into the Schrodinger equation Eq. (p^, multiplying 
by the operator e^'^^*-' from the Ihs, and using the fact that [^(t), S'(t')] = 
[S{t), S(t')] = 0, one obtains 

^^Mt)) + tS{t)mt)) = e-^W^(t)e^W|$(t)), (27) 

where the transformed Hamiltonian on the rhs can be expressed in terms of 
the commutator series (Baker-Campbell-Hausdorff expansion) 

e-^W^(t)e^W = H{t) + [H{t), S{t)] + ^[[H, S{t)], S{t)] + ■■■ (28) 

An important advantage of the CCE is that, due to the FS momentum re- 
strictions in Eq. (^61), the above series terminates after the first few terms 
(three for quartic interactions) and a closed-form expression of the trans- 
formed Hamiltonian ( pS]) can be obtained in terms of S{t). By requiring that 
all FS-pair creation processes are eliminated from the above equation, one 
obtains the CCE equations [0, |55| for S(t). Before proceeding with such 



a calculation however, one needs to derive explicit expressions for H{t) and 
U{t). 

4.3 Effective Hamiltonian and Transition Matrix Ele- 
ments 

4.3.1 Discussion for Doped Semiconductors 

In the general case, the effective Hamiltonian H{t) is given by Eq. (pi]) 
and the effective optical transition operator U(t) by Eq. (^). For calcu- 
lations, it is useful to re-express such equations in second-quantized form. 
In the case of undoped semiconductors, it was shown in Ref. [|154|| that 
the effective Hamiltonian H{t) describes the same interactions as the bare 
Hamiltonian H, however among quasiparticles with time-dependent proper- 
ties determined by the pump polarization. Similarly, the effective interband 
transition matrix element in U{t) was shown to have an imaginary part that 
describes the dephasing due to the exciton-exciton interactions. The purpose 
of this section is to present similar results in the case of a doped semicon- 



ductor 1 157, 156 
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We start with Eqs. (^Tj) and (P^]), which express H(t) and U(t) in terms 
of the canonical transformation operator 'P^(t). The latter is given by Eq. 
([23|), which includes the effects of the Coulomb interactions on the pump 
photoexcitation. It is important to realize that the effect of the e-h inter- 
actions on the pump and the probe photoexcitations is very different. For 
an adequate description of the FES, the e-h interactions should be taken 
into account non-perturbatively for the ]9rofe^photoexcited pair. Indeed, the 
nonlinear absorption spectrum at a given frequency u close to the FES reso- 
nance is determined by the time-evolution for long times (of the order of the 
dephasing time T2) of an e-h pair photoexcited by the probe at the energy 
u [|163|| . Since the characteristic "e-/i interaction time" E^l (inverse HFA 
bound state energy) that determines the non-exponential polarization decay 
of the FES is much shorter that the dephasing time, the long-time asymp- 
totics of the response function (to the probe) depends non-perturbatively 
on the e-h interactions. In contrast, a short pump optical pulse excites a 
wavepacket of continuum e-h pair states (unlike in the discrete exciton case) 
with energy width ~ t~^, which thus evolves during timescales comparable 
to the pulse duration tp. Also, the corrections to the effective Hamiltonian 
are determined by the time evolution of the pump-induced carriers only up 
to times ~ tp [see Eq. (pT|) ]. Therefore, if the "e-/i interaction time" E]^ is 
larger than tp, tpEM < 1, (i.e., if the pump pulse frequency width exceeds 
Em)i the e-h interactions can be treated perturbatively when describing the 
time-evolution of the pump-photoexcited pairs. This can also be shown ex- 
plicitly for the third-order nonlinear polarization. In the general expression 
for x^^\ ^ contributions that depend on the pump are integrated over the 
width of the pump pulse; therefore, any resonant enhancement of x^^^ that 
depends on the pump frequency will be broadened out for sufficiently short 
pulses with frequency width that exceeds Em- In other words, when deriv- 
ing the pump-renormalized parameters, one can treat Coulomb interactions 
perturbatively if the above condition is fulfilled. In fact, the above situa- 
tion is somewhat similar to the calculation of the linear absorption spectrum 
close to the indirect transition threshold, where perturbation theory can be 
used Thus the above consideration applies even for long pulse dura- 

tions provided that the detuning VL exceeds Em- However, in order to obtain 
the full absorption spectrum, the time-evolution of the ]5ro6e-photoexcited 
pair with such effective Hamiltonian (with perturbatively calculated time- 
dependent parameters) should be treated non-perturbatively. As can be seen 
from the above discussion, an important advantage of this formalism is that 
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it seprates naturally between the perturbative and the non-perturbative in- 
teraction effects. 



4.3.2 Second Quantization Expressions 

We now proceed with the derivation of the effective Hamiltonian. To lowest 
order in the interactions, V^{t) can be presented as 

q p,p',k 

+ ^ E n\(PP'; k; t)ai+p,_k&ip&lp'&-k, (29) 

P,p',k 

where Veh is the probability amplitude for excitation of an e-h pair with zero 
momentum satisfying 

i^Pe/.(q, t) = [n + £f ) + eL'i] VehicL, t) + fiSpit) - E ^(q - q')^e/.(q', t). 

q' 

(30) 

In Eq. (|2|), 

P|,(pp';k;t) = z f dt'e-^(*-*')(^+4+^p,-li+^;+p.-J 

X [Vip - k) [Pe/.(P', t') - Pe/.(P + P' - k, t')] - (P P')] (31) 

describes the scattering of the photoexcited e-h pair with an electron, and 
Pg\(pp';k;t) = -W dt'e V p+p'-"^ p' " V 

J — oo 

X [t;(p - k) [Vehip', t') - Vehip + P - k, t')] - (p ^ p')] (32) 
describes the scattering of the photoexcited e-h pair with a hole. The above 



expressions describe in the lowest order in the screened interaction [164 



?;(p — k) the coherent pump-induced processes, the effects of the Hartree- 
Fock pair-pair and pair-FS interactions, and the dynamical FS response to 
the pump photoexcitation. 

By substituting Eq. ( p9D into Eq. p^), we obtain the following expression 
for the effective optical transition operator: 

UKt)\Mt))= E Mp(t)at 61^10) 

p>kp 

+ \ E ^PP'k(t)«J.«>Lp-p'«k|0), (33) 

p,p'>kp,k<kp 
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where the effective matrix element Mp{t) includes corrections due to phase 
space filling and Hartree-Fock interactions, and Mpp/k(t) is the probability 
amplitude for indirect optical transitions induced by the pump optical 
field, which contribute to the pump-probe polarization in the second order in 
the interactions. The explicit expressions for Mp{t) and Mpp/k(t) are given 
in Appendix C. 

We turn now to the effective Hamiltonian H{t). After substituting Eq. 
(H) into Eq. (|l]) we obtain that 

(34) 



where 



= ^q ~ /^^p(^)Re 
is the effective conduction electron energy; 



Pe/^(q,t)-E^e/.(qq';q;^) 

q' 



is the effective valence hole energy; 

Veh{t) = 



kk'q 



is the effective e-h interaction; and 

^ kk'q 



(35) 



(36) 



(37) 



(38) 



is the effective e-e interaction. The explicit expressions for t;e/i(q; kk'; t) and 
t'ee(q; kk'; t) are given in Appendix C. As can be seen, liii) has the same 
operator form as the bare Hamiltonian H . However, both the effective band 
dispersions and the effective interaction potentials are now dependent on 
time. Note here that the above pump-induced renormalizations only last for 
the pulse duration tp. As discussed above, they are therefore perturbative in 
the screened interactions for tpEu < 1 or for > Em- 

Let us first discuss the effect of the pump-induced self-energy corrections 
to the conduction and valence band energies, given by the last terms in 
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Eqs. (^) and (PBI). The dispersion of the effective band is shown in Fig. 
|l[ As can be seen, the pump pulse leads to a bandgap increase as well as a 
change in the momentum dependence (band dispersion) that last as long as 
the pump pulse. The magnitude of the bandgap increase is of the order of 
{fi£pY/fl (for off-resonant excitation) and {fi£p)Hp (for resonant excitation) 
and leads to, e.g., the ac-Stark blueshift. As we shall see, for pulse duration 
shorter than the dephasing time, it also leads to bleaching and gain right 
below the onset of absorption, analogous to the case of excitons or two- 
level systems. It should be emphasized that these are coherent effects that 
should not be confused with the incoherent bandgap redshift due to the e-e 
interactions among real photoexcited carriers [3^. In particular, the above 
bandgap renormalization is induced by the transverse EM-field of the laser, 
as compared to the usual bandgap renormalization due to a longitudinal EM- 
field, i.e.. Coulomb screening. The change in the band dispersion, whose 
relative magnitude is of the order of {fiSp/fl)'^ (for off-resonant excitation) 
or [fiSptp)"^ (for resonant excitation), can be viewed as an increase in the 
effective density of states and, to the first approximation, in the effective 
mass. This is important in doped semiconductors because, as we shall see 
later, it leads to an optically-induced time-dependent enhancement of the 
e-h interactions and scattering processes with the FS electrons. 

The effective Hamiltonian H{t) also includes pump-induced corrections 
in the effective interaction potentials, determined by the pair-pair and pair- 
FS interactions during the pump photoexcitation. By expanding Eqs. 



and ( |129| ) for carrier energies close to the Fermi surface using Eqs. (|3T|) and 
([3^), one can show that these corrections vanish at the Fermi surface; for the 
typical FS excitation energies Ae ~ Em that contribute to the FES, their or- 
der of magnitude is {fiSpAe/fl"^)^ (for off-resonant excitation) or {/iSp Ae t^)^ 
(for resonant excitation). Thus the corrections to the interaction potentials 
are suppressed for below-resonant excitation by a factor of {Em/^Y-, or for 
short pulses by a factor of {EMtpY, as compared to the self-energy correc- 
tions. Such a suppression is due to the Pauli blocking effect and the screening, 
which leads to the vanishing of the pump-induced corrections to the inter- 
action potentials at the Fermi surface. Similarly, the pump-induced indirect 
optical transition matrix elements Mpp/k(t) are suppressed by the same fac- 
tor as compared to the direct transition matrix element Mp{t) [first term 
in Eq. (0)], while they contribute to the pump-probe polarization only in 
the second order in the screened interactions. Therefore, in the doped case, 
the screened Coulomb interaction leads to subdominant parameter renormal- 
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izations to the effective Hamiltonian H{t) and transition operator U^{t) for 
sufficiently short pump pulses or for off-resonant excitation. In contrast, 
in undoped semiconductors, the exciton-exciton interactions lead to non- 
perturbative contributions discussed in Ref. ||154|. 



4.4 Description of the Electron— Hole Pair Dynamics 



In this section, we present the final formulae for the nonlinear pump-probe 
polarization of the interacting system by applying the CCE to the effective 
Hamiltonian H{t) in order to treat the dynamical FS response. The CCE 
equation Eq. (|27|) contains the operator S(t) described by a hierarchy of cou- 
pled equations for the amplitudes s, S2, ■ ■ ■, defined by Eq. (PB]). As discussed 
at length above, in the coherent regime of interest here, the e-e scattering 
effects are suppressed, and the nonlinear absorption spectrum is dominated 
by the e-h interactions. This allows us to use the dephasing time approxima- 

m 



m 



tion for treating the probe-induced e-e scattering processes 
which case the above hierarchy terminates after s. Importantly, the e-h cor- 
relations (dynamical FS response) are still treated non-perturbatively, since 
they are determined by s |Q Then all FS-pair creation processes 

can be eliminated explicitely from the rhs of Eq. (^) , leading to the follow- 
ing nonlinear differential equation for the one-FS-pair scattering amplitude 
s(p,k,t): 



55, m, 



dt 



V 



.(t)-4(t)J s{p,k,t) 
+ ^(P',k,t) 

p'>kp 



k'<kp 



(39) 



The e-h scattering processes described by the above equation are sketched in 
Fig. Here V is the s-wave component ^ of the screened interaction 



approximated for simplicity by its value at the Fermi energy |^2|, ^ . 
This neglects plasmon effects, which are however small within the frequency 
range of the FES 
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Although this approximation is standard for 
the linear absorption case, its justification for the transient spectra requires 
more attention. Indeed, the characteristic time for screening buildup is of the 
where Qp is the typical plasma frequency corresponding to the FS 



order ^ , 



166| , |167| , |163|| . This time is however shorter than the typical pump duration 



lOOfs and dephasing time (which is of the order of ps near the Fermi 
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energy [^]), so the Coulomb interactions can be considered screened also 



for the nonlinear absorption case. In Eq. 
energy contribution to the excitation energy 



we neglected the hole recoil 
g since el^{t)-el{t) < Em 



due to a sufficiently heavy hole mass ||4^. Note that, by increasing the hole 
effective mass, the pump-induced hole self-energy, Eq. (|36D , reduces the hole 
recoil energy and thus the corresponding broadening ^2[. In real samples, 
the relaxation of the momentum conservation condition due to the disorder 
will also suppress the hole recoil broadening effects. In fact, the disorder can 
even lead to localized hole states |^5|], which corresponds to the infinite hole 
mass limit. 



From Eq. (27) we then easily obtain the following expression for 



p>kp 



(40) 



where is the creation operator of the zero-momentum hole state and the e- 
h pair wavefunction $p(t; t') satisfies the "Wannier-like" equation of motion, 

.d%{t,t') 



dt 



p'>kp 



where 



V{p,t) = V 



1 



(41) 



(42) 



k'Kkp 

is the effective e-h potential whose time- and momentum- dependence is de- 
termined by the response of the FS electrons to their interactions with the 
probe-induced e-h pair (vertex corrections) [sketched schematically in Fig. 
11(b), responsible for the unbinding of the HFA bound state]. 



ut) = v y: 



k'Kkp 



1+ E ^(P',k',t) 



p'>kf 



(43) 



is the self-energy due the to the sudden appearance of the photoexcited hole, 
which leads to non-exponential polarization decay [described by Imeyi(t)] 
due to the Anderson orthogonality catastrophe and a dynamical res- 
onance redshift [described by ReeAit)], and the dephasing width Fp de- 
scribes all additional dephasing processes (due to e-e interactions, hole re- 
coil, and phonons). Eq. (41) should be solved with the initial condition 
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$p(t',t') = Mp(t'), where the effective matrix element Mp{t') is defined by 
Eq. (0). 

It is worth stressing here the analogy between Eqs. ( ^OD and ( |41| ) and 
the corresponding problem in undoped semiconductors. Indeed, Eq. (^0]) is 
the direct analog of an exciton state, whereas again Eq. (^if) is very similar 



to a Wannier equation. However, Eqs. (|39D and (41) include the effects 
of the interactions between the probe-photoexcited e-h pair and the FS- 
excitations, and the wavefunction $p(t,t') describes the propagation of the 
photoexcited pair "dressed" by the FS excitations. Such a "dressing" is due 
to the dynamical FS response, which leads to the dynamical screening of 
the effective e-h interaction Eq. (^). The time-dependence of the latter is 
determined by the FS scattering amplitude s(p, k, t) and is affected by the 
pump excitation as described by Eq. (pOl) . One can easily verify that by 
setting s(p, k, t) = in Eq. (^Tj), we recover the results of the Hartree-Fock 



(ladder diagram, static FS |^) approximation. If one neglects the nonlinear 



(quadratic) term in Eq. (|39|) , one recovers the three-body (Fadeev) equations 
P2| , |168|| . Note that the coupled equations for $p(t, t') and s(p, k, t), obtained 
by neglecting the multipair excitations in Eq. (^Bf), can be extended to include 
the hole recoil-induced corrections ||58|| . 

Using Eqs. (^5]), (p6D , and (j^, we now can express the pump-probe 



polarization Eq. ( P^Uj ) in terms of the e-h wavefunction $ and the effective 
transition matrix element Mp. Assuming, for simplicity, a delta-function 
probe pulse centered at time delay r, Ss{t) = SsS{t — r), we obtain a simple 
final expression: 

Pk.(t) = te{t - r)/i2£,e^''-'--'^-(*-^) Mpm^it, t). (44) 

p>kp 

Eq. (P^ ) expresses the pump-probe polarization in terms of two physically 
distinct contributions. First is the effective transition matrix element Mp(t), 
which includes the effects of pair-pair and pair-FS interactions and Phase 
space filling effects due to the pump-induced carriers present during the 
probe photoexcitation. Second is the wavefunction $p(t) of the e-h pair 
photoexcited by the probe, whose time dependence, determined by Eqs. (^Tj) 
and (^), describes the formation of the absorption resonance. Despite the 
formal similarities, there are two important differences between the doped 
and the undoped cases. First, in the doped case, the time evolution of the 
e-h wavefunction $p(t) is strongly affected by the interplay between the 
e-h correlations and the pump-induced transient changes in the bandgap 
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and band dispersion relations. As we shall see later, this can be viewed 
as an excitation-induced dephasing. Second, unlike in the undoped case 
154|| , the pump-induced corrections in the effective matrix element Mp{t) 



are perturbative in the screened interactions if the pump detuning or the 
pump frequency width exceed the Coulomb energy Em- 

4.5 The Case of Two Coupled Subbands 

In this section, the above results are extended to the multisubband case in 
order to investigate the ultrafast PP dynamics of the FES-exciton hybrid 
formed in asymmetric QW's with partially occupied subbands ||169| , |17CI| , 



|T7I], |I7g, |T73|, |T7§. In such one-sided MDQW structures, interband 



optical transitions from the valence band to several conduction subbands 
are allowed due to the finite overlap between the hole and electron envelope 
wave-functions. When the exciton Fano resonance from a higher empty 
conduction subband is nearly resonant (within a few meV) to the Fermi 
level, the FES coming from the lowest occupied subband is enhanced by over 
two orders of magnitude [^, |169| , |170|| . Such many-body effects on the linear 



spectrum have been described by using the simple two-subband Hamiltonian 
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H = Y1 eik«lk«ik + 51(4 + ^g)^-k^-k - ^ij(q)4p+q«jp&k-q&k, (45) 

jk k ij pkq 

where al^^ and ef^ are the creation operator and the energy of a conduction 
electron in the ith subband, folj^. and are those of a valence hole {Eg is the 
bandgap), and fjj(q) is the screened e-h interaction matrix with diagonal 
(off-diagonal) elements describing the intrasubband (intersubband) scatter- 
ing. Due to the screening, the interaction potential is short-ranged and 
can be replaced by its s-wave component 0, close to the Fermi surface, 
Vij{q) ~ Vij ||171| , |172| , Here we consider the case where only the first sub- 



band is occupied, but the Fermi level is close to the exciton level (with binding 
energy Eb) below the bottom of the second subband. For large values of the 
FES-exciton splitting A — Ep — Eb, where A is the subband separation, the 
linear absorption spectrum consists of two well separated peaks, the lower 
corresponding to the FES from subband 1, and the higher corresponding to 
the Fano resonance from the exciton of subband 2 broadened by its coupling 
to the continuum of states in subband 1. With decreasing A — Ep — Eb, 
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the FES and the exciton become hybridized due to the intersubband scat- 
tering arising from the Coulomb interaction. This results in the transfer of 
oscillator strength from the exciton to the FES and a strong enhancement 
of the absorption peak near the Fermi level due to the resonant scattering of 
the photoexcited electron by the exciton level [|169| , |17CI| , |171i |172| , |173 |. 



Since the PP signal is linear in the probe field, the essential physics can be 
captured by assuming a 5-function probe pulse, Srit) = E^-e^'^^'^ 5{t — r), and 
a Gaussian pump pulse. The PP polarization has the form {t > r) 

P{t) = z^,e-*--*5]/i,/i,(0|f/i(t)/C(t,r)f/j(r)|0), (46) 

ij 

where /C(t, r) is the time-evolution operator for the effective Hamiltonian, 

Hit) = E 4k(i)«Uk + E 4W&lk&-k + v^hit) + vut), (47) 

jjk k 

where Veh and Vee are the effective e-h and e-e interactions and Uj{t) is the 
effective transition operator given below. Here 



4k(^) = + ^4k(^). = ^1 + ^ + A6i:(t) (48) 

are the band dispersions with pump-induced self-energies (to lowest order 
in the interactions) 

A4k(^) = -^p(t)Mk(i)+ Wk(t)]/2, 
Ael{t)= -£p{t)ReJ2^^^P^^it), (49) 



with Pik{t) satisfying 



+ el + fi)p,k(t) - E %P.q(^) + (50) 



9^ 

Note that the pump induces additional intersubband scattering, described by 
A^i2k(^)- To lowest order in the interactions, the effective transition operator 
appearing in Eq. (^) is given by 

t>/(t) = E</'^.k(t)a]k&Lk, (51) 

Jk 
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with 



1 



(52) 



takes a familiar 
in a 



In the one-subband case and in the coherent hmit, Eq 
form, 0k(^) = 1 — bk(^)P — the usual Pauli blocking factor ||157 
multi-subband case, the latter is a matrix. 

Similar to the one-subband case, the state /C(t, r)[//(r)|0), entering into 
(|46|), can be viewed as describing the propagation of the e-h pair (with 
wavefunction $jj(k, t)) excited by the probe pulse at time r, dressed by the 
scattering of the FS excitations (dynamical FS response). The latter leads to 
a dynamical broadening described by the amplitude Sij{p,k,t) that satisfies 
the differential equation [0, |157| 



dsij{p,k,t) 
dt 



+ E[^4p Wso(P, k, t) - Ae'i^^it)suip, k, t)] 
I 

I q>kp 



(53) 



with initial condition Sjj(p, k, r) = 0, and p and k labeling respectively the 
(ith subband) FS electron and the (jth subband) FS hole. Since only the 
first subband is occupied, the only non-zero components of Sij are Sii(p, k, t) 
and S2i(p, k, t), which describe the intra and intersubband FS excitations 
respectively. The photoexcited e-h pair wavefunction $jj(k, t, r) satisfies the 
Wannier-like equation 



' dt 



I 

- E Vii(k,t)<i>ij{q,t) 

l,q>kp 



jk(r), where 



(54) 



with initial condition $jj (k, r) = 

eA(t) = - E E sii(p',k',t)i;ii] (55) 

k'<kp p'>kp 

is the self-energy due to the readjustment of the FS to the photoexcitation 



of a hole |p.57| , pa^ and F is the inverse dephasing time due to all the processes 
not included in H. In Eqs. (53) and (^41), 
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Vij{k,t) = Vij - Sii{k,k',t)vij (56) 

l,k'<kp 

is the effective e-h potential whose time-dependence is due to the dynamical 
FS response |p.57|| . Note that it is the interplay between this effective potential 



and the pump-induced self-energies that gives rise to the unadiabatic FS 
response to the pump field. In terms of $jj(k, t), the polarization (^61) takes 
the simple form {t > r) 

Pit) = -^£,e-*-^*^/i,/i, Y: MKMjuit): (57) 

ijl k>kp 



with (/)jj(k, t) given by (|5^). The nonlinear absorption spectrum is then 
proportional to lmP{uj), where P{uj) is the Fourier transform of the rhs of 

5 Coherent Ultrafast Dynamics of the Fermi 
Edge Singularity 

5.1 Discussion of the Physics: Monochromatic Excita- 
tion 

In this section, the role of the e-h correlations in the nonlinear optical re- 
sponse of the FES is illustrated in a simple way for monochromatic excitation 
156| , [83| , |155| , |157|| . In the latter case, the theory discussed in the previous 



section applies for pump detunings larger than the characteristic Coulomb 
energy, Q > Em- As discussed above, close to the Fermi edge, the linear 
absorption spectrum of the FES can be approximately described using the 
simple power law expression in Eq. (^). The monochromatic pump excitation 
leads to a resonance blueshift, originating from the shift in the effective band 
energies [see Fig. 0], and to a bleaching mainly due to the Pauli blocking 
(Phase Space Filling) which reduces the effective transition matrix element 
(analogous to the dressed atom picture ||160| ). More importantly, however. 



the pump-induced change in the band dispersion increases the density of 
states close to the Fermi surface and thus also increases both the e-h scat- 
tering strength g = VM and the phaseshift 5 ~ tan~^(7r(7) that determine 
the FES lineshape (see Eq. (H)). This, in turn, leads to an increase in the 
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FES exponent (3 that determines the resonance width. In contrast, in the 
case of a bound excitonic state of dimensionahty D and Bohr radius a^, the 
resonance width remains unchanged, while the oscillator strength, oc a^^ , 
increases by a factor ~ (1 + D Am), where Am is the pump-induced change 
in the effective mass [|155| , p.56| , |83|| . Such an optically-induced enhancement 



of the exciton strength competes with the bleaching due to the Pauli blocking 
and the exciton-exciton interactions. This results in an almost rigid exci- 
ton blueshift, consistent with experiment ||175| , |176|| and previous theoretical 



results |T73, |178 



However, in the case of a FES resonance, the pump-induced change in 
the exponent 5 leads to a stronger oscillator strength enhancement than for 
a bound exciton state. Obviously, such an enhancement cannot be described 
perturbatively, i.e., with an expansion in terms of the optical field, since, as 
can be seen from Eq. (^), the corresponding corrections to the absorption 
spectrum diverge logarithmically for frequencies uj at the Fermi edge. Eq. 
d^) shows that the effect of the pump on the FES can be thought of as an 
excitation-induced dephasing that affects the frequency dependence of the 
resonance; again, this is in contrast to the case of the exciton. In the time 
domain, this also implies a memory structure related to the response-time of 
the FS-excitations. Therefore, the qualitative differences between the nonlin- 
ear optical response of the FES and the exciton originate from the fact that an 
exciton is a discrete bound state, while the FES is a continuum many-body 
resonance. The FS responds unadiahatically to the pump-induced change in 
the density of states via an increase in the e-h scattering of low-energy pair 
excitations. Such scattering processes, which determine the response of the 
Fermi sea to the hole potential in the course of the optical excitation, are 



responsible for the unbinding and broadening of the HFA bound state ||163 |. 
Therefore, the pump field changes the broadening and dephasing effects even 
for below-resonant photoexcitation. On the other hand, due to the finite 
Coulomb binding energy of the exciton, the pump optical field can polarize 
such a bound state and change its Bohr radius without ionizing it. 

5.2 Results for Short— pulse Excitation 

In this section, we discuss the nonlinear absorption of the FES in the case of 
short pulse excitation. The results presented here were obtained by solving 
numerically the differential equations (^T|) and (p9|), using the Runge-Kutta 



method, for Gaussian pulses with duration tp = 2.0Ep^ ||157|| . In order 
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to suppress the incoherent effects due to the e-e scattering of real carrier 
populations with the FS, below-resonant pump excitation and negative time 
delays were considered. Under such excitation conditions, the co/ierent effects 
in which we are interested dominate, and the Coulomb-induced corrections 
to the effective parameters, discussed in section are perturbative. The 
goal of this section is to study the role of the dynamical FS response [e-h 
correlations) on the pump-probe dynamics. For this reason, we compare the 
results of the theory outlined in section ^ to those of the HFA, obtained by 
setting s(p, k,t) = in Eq. (0). As mentioned above, in the latter case, the 
(spurious) HFA bound state does not interact with the FS pair excitations, 
even though it can merge with the continuum when one introduces a very 
short dephasing time. 

In Fig. ^, the linear absorption lineshape (i.e., in the absence of pump) of 
the FES is compared to the HFA (without the dynamical FS response). The 
parameter values g = 0.4 and F = 0.1£^i? were used, which were previously 
used to fit the experimental spectra in modulation doped quantum wells 
|]50| , For better visibility, we shifted the curves in order to compare 

their lineshapes. The linear absorption FES lineshape is consistent with that 
obtained in Ref. ||5^. On the other hand, the HFA spectrum is characterized 
by the coexistence of the bound state and a continuum contribution due 
to the fact that, in 2D, a bound state exists even for an arbitrary weak 
attractive potential. We note that if one limits oneself to linear absorption, 
it is possible to artificially shorten the dephasing time T2 = F"^, mainly 
determined by the hole recoil effects, by taking F ~ Em- Then the spurious 
discrete state and the continuum merge, and the discrepancy between the 
two linear absorption lineshapes decreases. This trick has been used for 
phenomenological fits of linear absorption experimental data. However, as 
we discuss below, in the nonlinear absorption case the differences in the 
transient spectra are significant so that the processes beyond HFA can be 
observed experimentally. 

Let us turn to the time evolution of the pump-probe spectra. In Fig. ^ we 
show the nonlinear absorption spectra calculated by including the dynamical 
FS response [Fig. ^(a)] and within the HFA [Fig. ^(b)] at a short time delay 
r = —tp/2. The main features of the spectrum are a pump-induced resonance 
bleaching, blueshift, and gain right below the onset of absorption. For off- 
resonant pump, these transient effects vanish for positive time delays after 
the pump is gone, and persist for negative time delays shorter than the 
dephasing time T2 = F"^. Similar features were also obtained for different 
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values of the pump amplitude, duration, and detuning. They are mainly due 
to the broadening induced by the transient renormalization of the energy 
band dispersion [Eqs. ( pSf ) and (|36|)] when its duration ~ tp is shorter than 
the dephasing time (analogous to excitons and two-level systems). 

Let us now turn to the role of the e-h correlations. In Fig. ^ we com- 
pare the differential transmission spectrum calculated by including the dy- 
namical FS response or within the HFA for long and short negative time 
delays. Note that, in PP spectroscopy, the experimentally measured dif- 
ferential transmission is given by Eq. (§) and in the weak signal regime, it 
reproduces the pump-induced changes in the probe absorption coefficient 
a{uj,T): DST{ijj,t) oc — Aa(c<j,r). Fig. ^(a) shows the results obtained for 
a long time delay, r = — I.5T2 = — 15.0i?^^, in which case frequency domain 
oscillations are observed. These oscillations are similar to those seen in un- 
doped semiconductors and two-level systems ||179||; however, their amplitude 



in the FES case is reduced. On the other hand, as shown in Fig. ^(b), for time 
delays comparable to the pulse duration, r = — O.IT2 = — tp/2 = — l.Oi?^^, 
the main features are a blueshift and bleaching. In this case the e-h corre- 
lations lead to a substantially larger width and asymmetric lineshape of the 
differential transmission spectrum. This comes from the different response 
of the FES to the pump-induced dispersion renormalizations when the e-h 
correlations are accounted for. This is more clearly seen in Fig. ^, where 
the magnitude of the resonance decrease, evaluated at the peak frequency, 
is plotted as a function of r. Clearly, the bleaching of the FES peak is sub- 
stantially stronger when the dynamical FS response is included than in the 
HFA case. Note that for |t| ~ F^^ the FES resonance is actually enhanced 
by the pump, as can be seen more clearly in Fig. |^ The time dependence of 
the resonance bleaching is strikingly different in the two cases. In the HFA 
case, the \DST{u,t) \ evaluated at the instantaneous peak frequency decays 
over a time scale |r| ~ F"^, i.e. during the dephasing time. This is similar to 
results obtained for a two-level system with the same effective parameters. 
On the other hand, the decay of \DST{(jj, r)| at the peak frequency is much 
faster when we take into account the e-h correlations. Note that the above 
results were obtained for off-resonant excitation. Under resonant conditions, 
a spectral hole is produced. In Fig. |^ we compare the resonance blueshifts, 
evaluated at the peak frequency, as a function of r. Again, a larger blueshift 
is predicted when the dynamical FS response is included. This suggests that 
in the experiment of Ref . |Q , where similar blueshifts were observed in two 
quantum well samples (one MDQW with a FES and one undoped sample 
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with a 2D-exciton) the effective parameters were larger in the latter case, 
due to the absence of screening and exciton-exciton interaction effects. 

In order to gain qualitative understanding of the role of the e-h cor- 
relations, let us for a moment neglect the momentum dependence of the 
pump-renormalization of the band dispersion and the phase space filling ef- 
fects and consider the bleaching caused by a rigid semiconductor band shift 
AEg(t), obtained from the pump-induced self-energies, Eqs. (|35D and (^61), 
evaluated at the bottom of the band (note that AEg{t) lasts for the duration 
of pulse). Within this model, the pump excitation has no effect on the e-h 
scattering amplitude s(p, k, t) [see Eq. (pQ])]. It is thus convenient to factorize 
the effects of the rigid band shift on the e-h wavefunction $p(t,t'): 

$p(t, t') = e-'It' '^^^'(*")'^*" |.p(t, t') . (58) 

This relation is general and defines $p(t,t'), which does not depend on 
AEgit). In the special case of a rigid shift, $p(t, t') coincides with $p(t — t') 
describing the propagation of the probe-photoexcited e-h pair in the absence 
of the pump pulse. By substituting into Eq. ( |58D the long-time asymptotic 
expression $p(t,t') = $p(t - t') oc [i{t - t')EFf~^ that gives the linear 
absorption spectrum of the FES at ^ Ep, and substituting the resulting 
$p(t,t') into Eq. (PD, we obtain a simple analytic expression for the effect 
of a pump-induced rigid band shift on the nonlinear absorption spectrum: 

aiu) oc Re dte*{-+*r)(t-r)-./; Ai?,(t')^t'[,(i _ (59) 

For AEgit) = one, of course, recovers the linear FES absorption in the 
vicinity of the Fermi edge ||163|| . For P = 0, Eq. ( p9|) gives the absorption of 
the non-interacting continuum. 

The physics of the FES can be seen from Eq. (|5^) . For /? = 1, this gives 
a discrete Lorentzian peak corresponding to the HFA bound state. However, 
during the optical transition, the e-h pair interacts with the FS electrons, 
leading to the readjustment of the FS density profile via the scattering of 
FS pairs. This results in the broadening of the discrete HFA bound state, 
which is governed by the time evolution of the FS. Such time evolution is 
unadiabatic due to the low-energy FS pairs, which leads to the character- 
istic power-law time dependence of the broadening factor in Eq. (|59|) . The 
interaction with the FS-pairs determines the exponent, < /3 < 1, of the 
latter, which leads to a non-Lorentzian lineshape in the frequency domain 
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and a non-exponential decay in the time domain. A detailed discussion of 
the above physics and the analogy to phonon sidebands and collision broad- 
ening may be found in Ref. ||163| ]. However, the CCE must be used in order 
to calculate the spectrum at all frequencies (and not just asymptotically 
close to Ep as with Eq. (|^)) and, most importantly, to describe the non- 
equilibrium FS and e-h pair response to the time-dependent increase in the 
effective mass/density of states, not included in Eq. (|59D . 

The resonance bleaching obtained from Eq. (^) as a function of r, is 
shown in Fig. ^(b) for j3 = 0.6, corresponding to the value of the param- 
eters used in Fig. |^(a) {g = 0.4), together with the HFA result (/3 = 1). 
Comparing Figs. ||(a) and (b), one can see that the rigid band shift approx- 
imation qualitatively accounts for the dynamics, but that there are strong 
discrepancies (see vertical scales), whose origin is discussed below. Both the 
magnitude and the time-dependence of the bleaching depends critically on 
the value of (3, which characterizes the interaction of the photoexcited e-h 
pair with the FS excitations. Because of such coupling, many polarization 
components are excited in the case of the continuum FES resonance, and it 
is their interference that governs the dynamics of the PP signal. Such in- 
terference is also responsible for the resonance enhancement and differential 
transmission oscillations at r < shown in Figs. |^ and ^. As (3 increases 
the interference effects are suppressed because the energy width of the con- 
tinuum states contributing to the FES narrows. In fact, this energy width 
is directly related to that of the linear absorption resonance. This is clearly 
seen in Fig. ^ where we show the effect of increasing g on the dynamics of 
the bleaching. It becomes more bound-state-like as, with increasing g, the 
FES resonance becomes narrower. On the other hand, in the HFA case, the 
decay rate is ~ T2 = F"^, i. e. it is independent on g, when Em becomes 
smaller than F, while for Em — F, the contribution of the continuum states 
produce a faster decay. 

Although the transient rigid band shift approximation, Eq. (|59D, explains 
some of the features of the dynamics of the bleaching, it strongly overesti- 
mates its magnitude. This is because Eq. (^) neglects the response of the 
many-body system to the pump-induced renormalization of the band's dis- 
persion. Such a transient change in the dispersion, which can be viewed as an 
increase in the density of states/effective mass for the duration of the pump, 
is important because it results in an enhancement of the e-h scattering. For 
example, in the case of monochromatic excitation, this leads to the change in 
the exponent (3 of the broadening prefactor in the integrand of Eq. (|59D , as 
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discussed Section Ol. For the short-pulse case, it is not possible to describe 



analytically the effect of the pump on the e-h scattering processes, due to 
the non-equilibrium unadiabatic FS response. The latter can however be 
described in a simple way with the numerical solution of Eqs. ( pQ]) and (^T]), 
which is presented in Figs. and |ll] and discussed below. 



In order to show the role of the pump-induced renormalization of the 
band dispersion in the presense of the dynamical FS response, we plot in 
Fig. |10| the function 

F{u,T) = \m J2 %i^,r) (60) 

p>kp 

where ^p{uj,t) is the Fourier transform of $p(t,r) defined by Eq. (^8|). 
Note that, in the presence of the band dispersion renormalization, the wave- 
function $p (which is independent of AEg{t)) no longer coincides with $p as 
in Eq. (|^. As can be seen in Fig. |Tn|(a), when the e-h correlations are taken 
into account, the pump-induced redistribution of oscillator strength between 
the states of the continuum that contribute to the resonance manifests itself 
as a dynamical redshift. This shift opposes the rigid band blueshift AEg{t) 
(when the latter is included). At the same time, the resonance strength is 
enhanced significantly. The latter effect originates from the interplay be- 
tween the transient increase in the effective mass/density of states of the 
photoexcited e-h pair and the "dressing" of this pair with the FS excitations 
[described by the effective potential V{p, t) in Eq. (|4l|) ]. In contrast, such an 
oscillator strength enhancement is suppressed in the HFA (which neglects the 
e-h correlations), as seen in Fig. |l3(b), in which case the main feature is the 
redshift of the resonance due to the pump-induced increase of the binding 
energy Em coming from the transient increase in the effective mass. 

In Fig. |TT] the effect of the renormalization of the band dispersion on 
the nonlinear absorption spectrum is shown. The optically-induced increase 
in the e-h interactions enhances significantly the strength of the FES and 
compensates part of the bleaching induced by the rigid band shift. A smaller 
enhancement is also seen in the HFA, where the pump-induced increase in the 
binding energy Em competes with the effects of the bandgap renormalization. 

5.3 Ultrafast Dynamics of the FES— Exciton Hybrid 

In this section we present the results for the evolution of the PP spectra 
of the FES-exciton hybrid [|180|| . The spectra in Ref. [|18CI|| were obtained 
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by the numerical solution of the coupled equations ( |5^ and (pH]), with the 
time-dependent band dispersions e^j^^{t) and e|^(t). In typical MDQW's, the 
intersubband scattering strength characterized by the screened potential vu 
is much smaller than the intrasubband scattering strength characterized by 
Vii (a value vu/vn ~ 0.2 was deduced from the fit to the linear absorption 
spectrum in ||171|| ). As discussed in the previous section, in the absence 
of coupling (t>i2 = 0), the different nature of the exciton and FES leads to 
distinct dynamics under ultrafast excitation. In the presence of the coupling, 
one should expect new effects coming from the interplay of this difference 
and the intersubband scattering that hybridizes the two resonances. Indeed, 
it was demonstrated in Ref. [ 18CI|| that, at negative time delays, the PP 
spectrum undergoes a drastic transformation due to a transient light-induced 
redistribution of the oscillator strength between the FES and the exciton. We 
discuss this at length in this section and show that such a redistribution is a 
result of the dynamical FS response to the pump pulse. In fact, the ultrafast 
PP spectra of the FES-exciton hybrid can serve as an experimental test of 
the difference between the FES and exciton dynamics. 



The calculations in Ref. ||18CI|| were performed at zero temperature for 



below-resonant pump with detuning f2 ~ Ep and duration tpEp/h = 2.0, 
and by adopting the typical values of parameters vu/vu = 0.2, P = OAEp, 
and viiAf = 0.3, Af being the density of states, previously extracted from 
fits to the linear absorption spectra ||171| , |43[| {Ep ~ 15 — 20 meV in typical 



GaAs/GaAlAs QW's [[7^, |S^). Note, however, that similar results were also 
obtained for a broad range of parameter values. In Fig. |12|(a) we plot the 
nonlinear absorption spectra at different negative time delays r < 0. For 
better visibility, the curves are shifted vertically with decreasing |r| (the 
highest curve represents the linear absorption spectrum). For the chosen 
value of the subband separation A, the FES and excitonic components of 
the hybrid are distinguishable in the linear absorption spectrum, with the 
FES peak carrying larger oscillator strength. It can be seen that, at short 
r < 0, the oscillator strength is first transferred to the exciton and then, 
with further increase in |r|, back to the FES. At the same time, both peaks 
experience a blueshift, which is larger for the FES than for the exciton peak 
because the ac-Stark effect for the exciton is weaker due to the subband 
separation A and the correlation effects. 

The transient exchange of oscillator strength originates from the different 
nature of the FES and exciton components of the hybrid. At negative time 
delays, the time-evolution of the exciton is governed by its dephasing time. 
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which is essentially determined by the homogeneous broadening F (in doped 
systems the exciton-exciton correlations do not play a significant role due to 
the screening). The pump pulse first leads to a bleaching of the exciton peak, 
which then recovers its strength at |r| ~ h/T. On the other hand, since the 
FES is a many-body continuum resonance, (i) the bleaching of the FES peak 
is stronger, and (ii) the polarization decay of the FES is determined not by 
F, but by the scattering with the low-lying FS excitations. This leads to 
much faster dynamics, roughly determined by the inverse Coulomb energy 
Em (see discussion in the previous section). However, the time-evolution 
of the hybrid spectrum is not a simple superposition of the dynamics of its 
components. Indeed, the pump-induced self-energies lead to the fiattening 
of the subbands or, to the first approximation, to a time-dependent increase 
in the effective mass (and hence the density of states), which in turn in- 
creases the e-h scattering | p.57|| . Important is, however, that, due to the 
subband separation and different nature of the resonances, such an increase 
is stronger for the FES. Therefore, the effect of the pump is to reduce the 
excitonic enhancement of the FES peak (coming from the resonant scattering 
of the photoexcited electron by the exciton level) as compared to the linear 
absorption case, resulting in the oscillator strength transfer from the FES 
back to exciton. In fact, such a transfer is strong even for smaller A [see Fig. 
|T2|(b)]. It should be emphasized that the above feature cannot be captured 
within the HFA. Indeed, the latter approximates the FES by a bound state 
and thus neglects the difference between the FES and exciton dynamics orig- 
inating from the unadiabatic response of the FS to the change in the e-h 
correlations. This is demonstrated in Fig. |l2|(c) where we show the spectra 
obtained without the FS dynamical response, i.e., by setting Sij = 0. Al- 
though in that case both peaks show blue shift and broadening, there is no 
significant transfer of oscillator strength 

This concludes our discussion of the coherent nonlinear response of the 
FES. In the rest of this article, we will review the role of size-dependent cor- 
relation effects, due to quasiparticle scattering via surface collective modes, 
on the ultrafast dynamics of the SP resonance in small metal nanoparticles. 
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6 Quasiparticle Scattering with Surface Col- 
lective Modes in Metal Nanoparticles 

6.1 Electron— Electron Interactions in Metal Nanopar- 
ticles 

Before we proceed with the ultrafast dynamics, we discuss in this section 
the effect of the surface collective excitations on the e-e interactions in a 
spherical metal particle. In particular, we present a detailed derivation of the 
dynamically screened Coulomb potential [p.52| , |153|| by generalizing a method 
previously developed for calculations of local field corrections to the optical 
fields jlSTl . 

The potential U{uj; r, r') at point r arising from an electron at point r' is 
determined by the equation 

U (to; r, r') = ^(r - r') + J dridr2u{r - ri)n(cj; ri, r2)U{uj; r2, r'), (61) 

where u{r — r') = e^|r — r'|^^ is the unscreened Coulomb potential and 
n(u;; ri, r2) is the polarization operator. There are three contributions to 11, 
arising from the polarization of the conduction electrons, the d-electrons, 
and the medium surrounding the nanoparticles: 11 = lie + lid + H^. It is 
useful to rewrite Eq. (pT|) in the "classical" form 



V- (E + 47rP) =47re25(r-r'), (62) 

where E{uj;r,r') = — Vf/(u;; r, r') is the screened Coulomb field and P = 
Pc + Pd + Pm is the electric polarization vector, related to the potential U 

as 

VP{uj; r, r') = -e^ / driU{uj; r, ri)f/(cu; ri, r'). (63) 

In the random phase approximation (RPA), the intraband polarization op- 
erator is given by 

ne(^;r,r') = E pe~f ^^i'!-n ^^(^)^g(^)^»*(^OC-(rO, (64) 



aa' a 



where and ip^ are the single-electron eigenenergies and eigenfunctions in 
the nanoparticle, and f{E) is the Fermi-Dirac distribution (we set h = 1). 
Since we are interested in frequencies much larger than the single-particle 



55 



level spacing, Ilc{uj) can be expanded in terms of l/u. For the real part, 
U'^^uj), we obtain in the leading order 

U'^iuj; r, n) = LvK(r) V5(r - n)], (65) 

where nc(r) is the conduction electron density. In the following we assume, 
for simplicity, a step density profile, nc(r) = nc9{R — r), where Uc is the 
average density. The leading contribution to the imaginary part, n'^'(co'), is 
proportional to cu"^, so that n^'(^) ^ ^'cil^)- 

By using Eqs. (^) and (|63D, we obtain a familiar expression for at 
high frequencies, 

P,(w; r, r') = ^^VU{uj; r, r') = e{R - r)xc{oo)E{uj; r, r'), (66) 

where Xc{^) = —e^fic/mu!^ is the conduction electron susceptibility. Note 
that, for a step density profile, Pc vanishes outside the particle. The d-hand 
and dielectric medium contributions to P are also given by similar relations, 

Pd{uj; r, r') = e{R - r)xd{ujMuj; r, r'), (67) 
P^uj; r, r') = e{r - R)Xm'E{iJ- r, r'), (68) 

where Xi = (ej~l)/4vr, i = d,m are the corresponding susceptibilities and the 
step functions account for the boundary conditions ||182|| . Using Eqs. (|66D- 
([68|), one can write a closed-form equation for U (cu; r, r'). Using Eq. (|63|) , the 
second term of Eq. (|6TD can be presented as — / (iriM(r — ri) V-P(co'; ri, r'). 
Substituting the above expressions for P, we then obtain after integration 
by parts 

|r — r I J |r — ri| 

+ eir-R)xm]U{u;ri,r') 



+ I [ dr^dr^-^ rn;'(cu; ri, r2)t/(cu; ra, r'), (69) 

J |r — ril 



with 

eicu) = 1 + Anxico) = ed{cu) - cul/cu^ (70) 

Up = Aire'^nc/m being the plasmon frequency in the conduction band. The 
last term in the rhs of Eq. (|69|) , proportional to Il"{uj), can be regarded as a 
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small correction. To solve Eq. we first eliminate the angular dependence 
by expanding U{uj;r,r') in spherical harmonics, Ylm{^), with coefficients 
Ulm{^'i r'). Using the corresponding expansion of |r— r'|~^ with coefficients 
QLHi'T^r') = 21^1'^"^"^'^'^ ^ ^ following equation for 



L -\- \ f T \ ^ 

e{u)ULM{^] r, r') = QiM^r, r') + 47r [x{uj) - Xm\ l^J—[ y^j Ulm{^] R, r') 

L'M' 

X UL'M'{uj;r2,r'), (71) 

where 

^LM,L'M'i^'^n,r2) = J rffi(ir2FiV(ri)n"(t^;ri,r2)yL'M'(r2), (72) 

are the coefficients of the multipole expansion of n^'(co'; ri, r2). For 11^' = 0, 
the solution of Eq. ( [7T| ) can be presented in the form 

ULMiuj; r, r') = a(^)e2Qi,,(r, r') + 6(^)^^-^, (73) 

with frequency-dependent coefficients a and b. Since n"(ci;) ^ n'^(ti;) for the 
relevant frequencies, the solution of Eq. (0) in the presence of the last term 
can be written in the same form as Eq. (|73D , but with modified a{uj) and 
b{uj). Substituting Eq. ( [73D into Eq. (|7ID, we obtain after lengthy algebra in 
the lowest order in U" 

aiu) = e-\u;), b{u;) = e^i(^) - e-\u;), (74) 

where 

is the effective dielectric function, whose zero, e'i{ujL) = 0, determines the 



frequency of the collective surface excitation with angular momentum L ||103 
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In Eq. (^), characterizes the damping of the L-pole collective mode 

by single-particle excitations, and is given by 

= ^ E \M!^^'?\f{K) - fiK'MK - K' + (77) 

where M^^ are the matrix elements of r^Y^Mi^)- Due to the momentum 
non-conservation in a nanoparticle, the matrix elements are finite, which 



leads to the size-dependent width of the L-pole mode |p.07|, |181 



7l = ^e'^d^). (78) 

L UJp 

For uj ~ ul, one can show that the width, 7^ ~ vp/R, is independent of 
Note that, in noble metal particles, there is an additional d-electron 
contribution to the imaginary part of 6^(0;) at frequencies above the onset A 
of the interband transitions. 

Putting everything together, we arrive at the following expression for the 
dynamically-screened interaction potential in a nanoparticle: 

, u{r — r') ^ 47r 1 




with e^^(co') = ej^^(co') — e^^{uj). Equation ([T9|) represents a generalization 
of the plasmon pole approximation to spherical metal particles. The two 
contributions to the rhs originate from two types of dynamical screening. 
The first describes the usual bulk-like screening of the Coulomb potential by 
the electrons inside the particle. The second contribution is a new effective 
interaction induced by the surface: the potential of an electron inside the 
nanoparticle excites high-frequency surface collective modes, which in turn 
act as image charges that interact with the second electron. It should be 
emphasized that, unlike in the case of the optical fields, the surface-induced 
dynamical screening of the Coulomb potential is size-dependent. 

Note that the excitation energies of the surface collective modes are lower 
than the bulk plasmon energy, also given by Eq. ( [76D but with = 0. As 
discussed below, this opens up new channels of quasiparticle scattering. 
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6.2 Conduction Electron Scattering 

In this section, we present the calculation of the rates of electron scattering in 
the conduction band accompanied by the emission of surface collective modes 
and discuss its possible experimental manifestations. In the first order in the 
surface-induced potential, given by the second term in the rhs of Eq. ([79D, 
the corresponding scattering rate can be obtained from the Matsubara self- 
energy 11] 



1 47re2 IM^^P 



where = {iuj — E^)~^ is the no n- interacting Green function of the con- 
duction electron. Here the matrix elements M^^f are calculated with the 
one-electron wave functions ^Jaij) = Rni{r)Yim{T) . Since |a) and \a') are 
the initial and final states of the scattered electron, the main contribution 
to the Lth term of the angular momentum sum in Eq. (|80|) will come from 
electron states with energy difference Ea — Ea' ~ ujl- Therefore, can 
be expanded in terms of the small parameter Eo/\E^ — E^,\ ~ Eq/ul, where 
Eq = {2mB?)~^ is the characteristic confinement energy. The leading term 



can be obtained by using the following procedure ||107| , |181|| . We present 
AC as 

M- = (c.a|r-yi„(r)|c.a'> = (c. 1^. Wr)ll|c.a-> _ 

where H = Hq + V{r) is the Hamiltonian of an electron in a nanoparti- 
cle with confining potential V{r) = Vo9{r — R). Since [H,r'"YLM{^)] = 
— ^V[r^yiA/(r)] ■ V, the numerator in Eq. (^Tj) contains a term proportional 
to the gradient of the confining potential, which peaks sharply at the surface. 
The corresponding contribution to the matrix element describes the surface 
scattering of an electron making the L-pole transition between the states 
|c, a) and |c, a'), and gives the dominant term of the expansion. Thus, in the 
leading order in \E^ — E^,\~^, we obtain 

j^LM ^ {c,aMr'^YLMir)]-\/Vir)\c,a') 
m{El - El,Y 



m{E-^-El^,Y 



VoRnl{R)Rn'l'{R)^lm.l'm'J (§2) 
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with ypf^,^, = JdvYilir)YLMim'n.'{^). Note that, for L = 1, Eq. 
becomes exact. For electron energies close to the Fermi level, E!^i ~ Ep, the 
radial quantum numbers are large, and the product VoRni{R)Rn'i'{R) can be 
evaluated by using semiclassical wave-functions. In the limit Vq —>■ oo, this 



product is given by flO^ 2JE^iE^,i,/R^, where E^^ = 7i'^{n + l/2fEo is the 



electron eigenenergy for large n. Substituting this expression into Eq. 
and then into Eq. (^OD, we obtain 



X ^^^^G^,K + zu;), (83) 



with 



Civ = E I^LVr = ^ if ^ f dxP,{x)P,{x)Mx), (84) 

M,m' 

where Pi{x) are Legendre polynomials; we used properties of the spherical 
harmonics in the derivation of Eq. (|8^ . For E'^^ ~ Ep, the typical angular 
momenta are large, / ~ fc^i? ^ 1, and one can use the large-/ asymptotics of 
Pi] for the low multipoles of interest, L <^ I, the integral in Eq. ([8^ ) can be 
approximated by 2F+i^ii'- After performing the Matsubara summation, we 
obtain for the imaginary part of the self-energy that determines the electron 
scattering rate 

where N{E) is the Bose distribution and gi{E) is the density of states of a 
conduction electron with angular momentum /, 



gi{E) = 2Y^6iE^^,-E)c^^^^-^, (86) 

where we replaced the sum over n by an integral (the factor of 2 accounts 
for spin). 

Each term in the sum in the rhs of Eq. (p5| ) represents a channel of electron 
scattering mediated by a collective surface mode with angular momentum L. 
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For low L, the difference between the energies of modes with successive values 
of L is larger than their widths, so the different channels are well separated. 
Note that since all ljl are smaller than the frequency of the bulk plasmon, 
one can replace 61,(0;) by 6^(0;) in the integrand of Eq. (|85|) for frequencies 

(J ~ UJl- 

Consider now the L = 1 term in Eq. (^), which describes the SP- 
mediated scattering channel. The main contribution to the integral comes 
from the SP pole in e5~^(a;) = 3e~^{uj), where es{uj) is the same as in Eq. (pIS)). 
To evaluate the integral in Eq. (|85|) , we can in the first approximation replace 
lmej^{uj) by a Lorentzian, 



Im 



(87) 



where Ug = uji = Up/ \J e'^iug) + and 7 = 7s + Ust'l^ius) are the SP 
frequency and width, respectively. For typical widths 7 <^ cUs, the integral 
in Eq. (|85|) can be easily evaluated, yielding 



Using the relation e^kple'^i^ioJs) +2em]~"^ = 'i'noj'll^EF, the SP-mediated scat- 
tering rate, ■y^{E^) = — ImE^(£^^), takes the form 

YAE) = 9.^4- ( [1 - fiE - -.)]. (89) 



uOg Ep \ E^ 

Recalling that Eq = {2mR'^)^^, we see that the scattering rate of a conduction 
electron is size-dependent 7I oc i?^^. At E = Ep + uJs, the scattering rate 
jumps to the value 97r(l + Us/Ef)Eq/uJs, and then increases with energy 
as E^^"^ (for cug <C Ep). This should be contrasted with the usual (bulk) 
plasmon-mediated scattering, originating from the first term in Eq. (|79D, in 
which case the rate decreases as E'^/"^ above the onset PS . 



Note that the total electron scattering rate is the sum, 7e + 7e, of the 
SP-mediated (7I) and the bulk-like (7e) scattering rates. In order to be 
observable, the former should exceed the latter. The typical size at which 
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7e becomes important can be estimated by equating 7I and the Fermi liquid 
e-e scattering rate [||], 7e(-E) = -j^k^^^-^f^- For energies E ~ Ep + Ug, the 
two rates become comparable for 

(kRf ^12^(1 + ^)"'(J^Y\ (90) 



In the case of a copper nanoparticle with uog — 2.2 eV, we obtain k^R ~ 8, 
which corresponds to a radius of i? ~ 3 nm. At the same time, in this energy 
range, the width 7I exceeds the mean level spacing 5, so that the energy 
spectrum is still continuous. The strong size dependence of 7I indicates 
that, although 7^ increases with energy slower than 7e, the SP-mediated 
scattering should dominate for nanometer-sized particles. Note also that the 
size and energy dependences of the scattering in the different channels are 
similar: the rate of scattering via the Lth channel is given by Eq. ( ^9]) with ujg 
replaced by dj^, Eq. ([76|), and the numerical factor 9 replaced by 3L(2L + 1). 

Concluding this section, we have shown that the SP-mediated scattering 
is the dominant scattering mechanism of conduction electrons in nanometer- 
sized nanoparticles for energies larger than ujs but smaller than ujp. The 
scattering rate in the Lth channel, 7^', increases with energy, in sharp con- 
trast with the bulk-plasmon-mediated scattering rate. The total scattering 
rate as a function of energy represents a series of steps aX E = uj^, on top of 
a smooth energy increase. We expect that this new effect should be observ- 
able experimentally by measuring e-e scattering rate in size-selected cluster 



beams in time-resolved two-photon photoemission spectrum [183 



6.3 Band Hole Scattering 

We now turn to the interband processes in noble metal particles and consider 
the scattering of a d-hole into the conduction band. The corresponding 
surface-induced potential, given by the Lth term in Eq. ([79D, has the form 

r.r') = 5^^^ (^)V.M(f)n-„(r'). (91) 

where 6^(0;) is given by Eq. ( [r5|) With this potential, the d-h.o\e self-energy 
is given by 
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where M^^ = {c,a\r^YLM{^)\d,<y') is the inter^anc? transition matrix ele- 
ment [compare with Eq. (0)]. Since the final state energies in the conduc- 
tion band are high (in the case of interest here, they are close to the Fermi 
level), the matrix element can be approximated by a bulk-like expression 
M^^f = 6aa'M^^^ , the corrections due to surface scattering being suppressed 
by a factor of {kpR)^^ ^ 1. 

The largest contribution to the self energy Eq. (|92D comes from the dipole 
channel, L = 1, mediated by the SP. In this case, after performing the 
frequency summation, we obtain for ImS^ 



where E"^ = E^ — E^ and /x is the interband dipole matrix element | ]152| | . We 
see that the scattering rate of a d-hole with energy E^, •y^E^) = ImS^(i?^), 
has a strong dependence on the nanoparticle size, which is, however, 
different from that of the intraband scattering, Eq. (|89|). 

The most important difference between the interband and the intraband 
SP-mediated scattering rates hes in their dependence on energy. Since the 
surface-induced potential, Eq. (^I]), only allows for vertical (dipole) inter- 
band single-particle excitations, the phase space available for the scattering 
of a (i-hole with energy E^ is restricted to a single final state in the con- 
duction band, with energy E^. As a result of this restriction, the (i-hole 
scattering rate, •yf^^E^), exhibits a peak as the difference between the energies 
of final and initial states, E^ = E^ — E^, approaches the SP frequency Ug 
[see Eq. (pS])]. In contrast, the energy dependence of 7| is smooth due the 
larger phase space available for scattering within the conduction band. This 
leads to the additional integral over final state energies in Eq. (|53|), which 
smears out the SP resonant enhancement of the intraband scattering. 

As we show later, the fact that the scattering rate of a d-hole is dominated 
by the SP resonance affects strongly the nonlinear optical dynamics in small 
nanoparticles. This is the case, in particular, when the SP frequency, Ug, 
is close to the onset of interband transitions. A, as, e.g., in copper and 
gold nanoparticles [|103|, |126|, |127||. Indeed, if the optical pulse excites an 



e-h pair with excitation energy u close to A, the d-hole can subsequently 
scatter into the conduction band by emitting a SP. According to Eq. (p3D, 
for uj ^ ujg, such a scattering process should be resonantly enhanced. In 
order to have an observable effect on the absorption spectrum, the scattering 
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rate of the photoexcited d-ho\e should be comparable (or larger) than that 
of the photoexcited electron. Close to Ep, the electron scattering in the 
conduction band comes from a two-quasiparticle process; the corresponding 

-2 eV. 



rate in noble-metals is estimated as ||183|| ■je ~ 10~ eV. If one assumes the 
bulk value for fi ~ 1 eV near the L-point ||184|| ), then 7^ exceeds 

7e for R < 2.5 nm. In fact, one would expect that, in nanoparticles, is 
larger than in the bulk due to the localization of the conduction electron 
wave-functions 11103 



7 Surface Plasmon Nonlinear Optical dynam- 
ics 

In this section, we discuss the effect of the SP-mediated interband scattering 
on the ultrafast optical dynamics in noble metal nanoparticles. We are inter- 
ested in the situation when the hot electron distribution has already thermal- 
ized and the electron gas is cooling to the lattice (stage III). In this case the 
transient response of a nanoparticle can be described by the time-dependent 
absorption coefficient a{u!,t), given by Eq. (|12]) with time-dependent tem- 



perature. In noble-metal particles, the temperature dependence of a origi- 
nates from two different sources. First is the phonon-induced correction to 
7s, which is proportional to the lattice temperature Ti{t). As mentioned in 
section ^, for small nanoparticles this effect is relatively weak. Note that, 
as discussed e.g. in Ref. [ [L37|] , the contribution to js coming from the e-e 
interaction ||125|, |129|, |130|, |127|, |139|, |14CI||, which depends on the electron tem- 



perature, plays a minor role when the interband transitions are resonantly 
excited. Second, near the onset of the interband transitions. A, the absorp- 
tion coefficient depends on the electron temperature T{t) via the interband 
dielectric function ed{uj) [see Eqs. (0) and (|13])]. In fact, in copper or gold 
nanoparticles, Ug can be tuned close to A, so the SP damping by interband e- 
h excitations leads to an additional broadening of the absorption peak ||103 



In this case, it is the temperature dependence of ed{uj) that dominates the 
pump-probe dynamics. Below we show that, near the SP resonance, both the 
temperature and frequency dependence of ed{uj) = 1 + iiiXdi^^) are strongly 
affected by the SP-mediated interband scattering. 

We start with the RPA expression for the interband susceptibility p3| , 
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XS^) = - E E Gt{iuj')G'^^{iu' + z^), (94) 

where G^(za;') is the Green function of a (i-electron. With the d-band fully 
occupied, the only allowed SP-mediated interband scattering is that of the 
d-hole. We assume here, for simplicity, a dispersionless d-band with energy 
E^. Substituting Gi{iuj') = [iuj' - E'^ + Ep - T.i{iuj')]-\ with T^Oiu) given 
by Eq. and performing the frequency summation, we obtain 

, eV r dE^gjE^) fjE^)-! 

^^^^ m2 J {E^dy uj-E'^ + iYh{uJ,E-y ^ ' 

where g{E^) is the density of states of conduction electrons. The scattering 
rate of a d-hole, 7j^(u;,i?'^) = ImS'^(£"^ — cj), is obtained from Eq. (|93D with 
Ed = E^- u: 



where Nioj) is negligible for frequencies ~ ^ /^b^" ||152|| . The rate 
7^(1^7, i?^) exhibits a sharp peak as a function of the frequency of the probe 
optical field. The reason for this is that the scattering rate of a d-hole 
with energy E depends explicitly on the difference between the final and 
initial states, E'^ — E, as discussed above; therefore, for a d-hole with energy 
E = E^ — uj, the dependence on the final state energy, E'^, cancels out in 
es{E'^—E) [see Eq. (0)]. In other words, the optically-excited d-ho\e scatters 
resonantly into the conduction band as u approaches Ug. It is important to 
note that 7^(0;, E'^) is, in fact, proportional to the absorption coefficient a{u) 
[see Eq. (0)]. Therefore, a and 7^ should be calculated self-consistently from 
Eqs. (0), (0), (|9|), and (H). 

It should be emphasized that the effect of 7^ on e'^i^uj) increases with 
temperature. Indeed, it can be seen from Eq. (p6| ) that the value of 7^ is 
appreciable only if E'^ — Ep > ksT. Since the main contribution to x'^ioo) 
comes from energies i?'^ ~ u; — A + Ep, the effect of ci-hole scattering on 
the absorption becomes important only for elevated electron temperatures: 
ksT > Ug — A. As a result, near the SP resonance, the time evolution of the 
differential absorption, which is governed by the temperature dependence of 
a, becomes strongly size-dependent. 
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The numerical results discussed below are taken from Refs. |T^, 153|| . In 



the experiment of Bigot et. al. ||126|| , the pump-probe measurements were 



performed on i? ~ 2.5 nm copper nanoparticles. The SP frequency, ujg — 2.22 
eV, was slightly above the onset of the interband transitions, A ~ 2.18 eV. In 
order to describe the time-evolution of the differential absorption spectra, one 
first needs to determine the time-dependence of the electron temperature, 
T(t), due to the relaxation of the electron gas to the lattice. For this, a 
simple two-temperature model is employed, defined by heat equations for 
T{t) and the lattice temperature T/(t): 

C(r)^ = -G{T-T,\ 

= G{T-Ti), (97) 

where C{T) = 7T and Ci are the electron and lattice heat capacities, re- 
spectively, and G is the e-p coupling [|141| , |142|, |143|| . The parameter values 



used in Refs. |T5|, |5|] were G = 3.5 x 10^^ Wm~^K-\ 7 = 70 Jm^^j^-^, 
and Ci = 3.5 Jm~^K~^, and the initial condition was taken as Tq = 1000 
K. The time-dependent absorption coefficient a{u!,t) was calculated self- 
consistently; the differential transmission is proportional to ar{uj) — a{u!,t), 
where ar{uj) was calculated at the room temperature. 

In Fig. |l^, the calculated nonlinear absorption spectra are shown for 
various nanoparticle sizes. Fig. |T3|(a) shows the spectra at several time delays 
for R = 5.0 nm; for this size, the SP-mediated c?-hole scattering has no 
effect. With decreasing nanoparticle size, the linear absorption spectra are 
not significantly altered, as can be seen in Figs. |13|(b) and (c). However, 
the change in the nonlinear absorption spectra becomes pronounced at short 
time delays corresponding to higher temperatures [see Figs. |r^(b) and (c)]. 
This effect is more clearly seen in the differential transmission spectra, shown 



in Fig. |T^, which undergo a qualitative transformation with decreasing size. 

Note that it is necessary to include the intraband e-e scattering in order to 
reproduce the differential transmission lineshape observed in the experiment 



126 |. For optically excited electron energy close to Ep, this can be achieved 
by adding the e-e scattering rate ||] -feiE") oc [1 - f{E'')][{E^ - Epf + 
(tt/cbT)^] to 7^ in Eq. (pSj). The difference in 7e(-E^) for E'^ below and above 
Ep leads to a lineshape similar to that expected from the combination of 
red-shift and broadening [see Fig. |l^(a)]. 

In Figs. p!4|(b) and (c) the differential transmission spectra are shown 
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with decreasing nanoparticle size. For R = 2.5 nm, the apparent red-shift 
is reduced [see Fig. |l^(b)]. This change can be explained as follows. Since 
here ~ A, the SP is damped by the interband excitations. This broadens 
the spectra for a; > a;^, so that the absorption peak is asymmetric. The d- 
hole scattering with the SP enhances the damping; since the w-dependence 
of 7^ follows that of a, this effect is larger above the resonance. On the 
other hand, the efficiency of the scattering increases with temperature, as 
discussed above. Therefore, for short time delays, the relative increase in the 
absorption is larger for u > Us- With decreasing size, the strength of this 
effect increases further, leading to an apparent blue-shift [see Fig. |l^(c)]. 
Such a strong change in the absorption dynamics originates from the 
dependence of the d-ho\e scattering rate; reducing the size by the factor of 
two results in an enhancement of 7^ by an order of magnitude. 

In Fig. |T3|, the time evolution of the differential transmission are shown 
for several frequencies close to Ug- It can be seen that the relaxation is slow- 
est at the SP resonance; this characterizes the robustness of the collective 
mode, which determines the peak position, versus the single-particle excita- 
tions, which determine the resonance width. For larger sizes, at which 7^ is 
small, the change in the differential transmission decay rate with frequency 
is smoother above the resonance [see Fig. |l^(a)]. This stems from the asym- 
metric lineshape of the absorption peak, mentioned above: the absorption is 
larger for uj > Ug, so that its relative change with temperature is weaker. For 
smaller nanoparticle size, the decay rates become similar above and below Us 
[see Fig. |15|(b)]. This change in the frequency dependence is related to the 
stronger SP damping for u > Us due to the d-ho\e scattering, as discussed 
above. Since this additional damping is reduced with decreasing tempera- 
ture, the relaxation is faster above the resonance. This rather "nonlinear" 
relation between the time-evolution of the pump-probe signal and that of 
the temperature becomes even stronger for smaller sizes [see Fig. |l5|(c)]. In 
this case, the frequency dependence of the differential transmission decay be- 
low and above Ug is reversed. Note that a frequency dependence consistent 
with our calculations presented in Fig. |TS|(b) was, in fact, observed in the 



experiment of Ref. [ 126 |, shown in Fig. |T5|(d). 
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8 Conclusions 



In this article we presented an overview of ultrafast nonlinear optics exper- 
iments that probe the dynamics of modulation-doped quantum wells and 
metal nanoparticles. We also discussed some examples where the ultrafast 
dynamics of such confined Fermi seas cannot be described in terms of the 
usual dephasing and relaxation time approximations. These examples show 
the important role of dynamical and non-equilibrium many-body correla- 
tions during femtosecond time scales. 

More specifically, we discussed a recent theory for the ultrafast nonlinear 
optical response of the FES. We focusscd on coherent effects, which dominate 
the pump-probe spectra during negative time delays and off-resonant exci- 
tation conditions. We demonstrated that the dynamical FS response leads 
to quahtatively different coherent dynamics of the FES as compared to the 
Hartree-Fock approximation. In the latter case, the time evolution of the 
resonance bleaching is governed by the dephasing time. In contrast, in the 
former case, polarization interference and e-h correlation effects dominate. 
This results in an initial fast FES dynamics, with a response time determined 
by the characteristic inverse Coulomb energy Em as well as an apparent res- 
onance enhancement during negative time delays, followed by a long-time 
decay determined by the dephasing time. Such dynamical features should be 
observable in ultrafast PP experiments. Using a simple model, we showed 
that the different dynamics of the FES and Hartree-Fock treatment can be 
attributed to the non-Lorentzian broadening of the HFA bound state due to 
its interactions with the gapless FS excitations, a process which is, of course, 
beyond the dephasing time approximation. In addition, we showed that the 
pump excitation directly affects the strength of the e-h scattering processes, 
which changes the frequency dependence of the resonance. The latter can be 
thought of as an excitation-induced dephasing effect that leads to a transient 
enhancement of the FES. These results indicate that ultrafast spectroscopy 
provides a powerful tool to study the role of correlations in the nonlinear re- 
sponse of a Fermi liquid during time scales shorter than the dephasing times. 
The dynamical features discussed above can also be used as an experimental 
signal to probe the crossover from FES to exciton bound state (exciton Mott 
transition) as a function of the FS density. 

The above e~h correlation effects and the fundamental differences in the 
dynamics between an exciton state and a FES many-body resonance can 
be best observed experimentally in asymmetric modulation-doped quantum 
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wells in the case where the lowest occupied subband is in close proximity to 
the second unoccupied subband and is coupled strongly to it via intersub- 
band e-h scattering. In such a system there is a strong redistribution of the 
oscillator strength between the FES and exciton peaks which is caused by 
the different dynamics of the FES and exciton components of the hybrid as 
well as by their coupling due to the e-h correlations. This originates from 
the dynamical Fermi sea response and leads to a strong transient changes in 
the PP spectra. 

We also discussed the role of size-dependent correlations in the electron 
relaxation in small metal particles. We identified a new mechanism of quasi- 
particle scattering, mediated by collective surface excitations, which origi- 
nates from the surface-induced dynamical screening of the e-e interactions. 
The behavior of the corresponding scattering rates with varying energy and 
temperature differs substantially from that in the bulk metal. We showed 
that the conduction electron scattering rate increases with energy, in sharp 
contrast to the bulk behavior, which could be observed in two-photon pho- 
toemission measurements. We also found that in noble metal particles, the 
resonant energy dependence of the d-hole scattering rate affects strongly the 
differential absorption. 

An important aspect of the SP-mediated scattering is its strong depen- 
dence on size. Our estimates show that it becomes comparable to the usual 
Fermi-liquid scattering in nanometer-sized particles. This size regime is, 
in fact, intermediate between "classical" particles with sizes larger than 10 
nm, where the bulk-like behavior dominates, and very small clusters with 
only dozens of atoms, where the metallic properties are completely lost. 
Although the static properties of nanometer-sized particles are also size- 
dependent, the deviations from their bulk values do not change the qual- 
itative features of the electron dynamics. In contrast, the size-dependent 
many-body effects discussed here do affect the dynamics in a significant way 
during time scales comparable to the relaxation times. The SP-mediated 
interband scattering reveals itself in the size-dependence of the transient 
pump-probe spectra. In particular, as the nanoparticle size decreases, the 
calculated time-resolved differential absorption lineshape shows a transition 
from an apparent redshift to a blueshift. This transition, absent in the RPA, 
comes from the correlations between collective surface and single-particle 
excitations. At the same time, near the SP resonance, these correlation leads 
to significant size-dependent changes in the frequency dependence of the re- 
laxation time of the pump-probe signal. These results indicate the need for 
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a systematic experimental studies of the size-dependence of the transient 
nonhnear optical response, as we approach the transition from boundary- 
constrained nanoparticles to molecular clusters. We expect that in the com- 
ing years ultrafast nonlinear optical spectroscopy will provide new insight 
into many-body effects in strongly correlated, and magnetic Fermi sea sys- 
tems psi wm, m M itm wm. 
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A APPENDIX A 



In this appendix we briefly outline the formalism of Ref . ||154|| . The pump- 
probe signal is determined by the polarization, 

P{t) = ^e-*'^^*(^(t)|f/|^(t)), (98) 

where the state |^E'(t)) satisfies the time-dependent Schrodinger equation. 



.d_ 



H - Hp{t) - Hs{t) 



\^{t)) = 0, (99) 



with the Hamiltonians H and Hp^s{t), given by Eqs. (|T6|)-(p!7D. The Hilbert 
space of the bare semiconductor, i.e. in the absence of optical fields, consists 
of disconnected subspaces C{^eh} which are labeled by the number of (in- 
terband) e-h pairs, I'eh- The corresponding bare Hamiltonian, H, conserves 
the number of e-h pairs in each band separately and in the Cl'^e/tj-basis 
has a block-diagonal form. The Hamiltonians Hp^s{t) couple the different 
subspaces Ci^eh} by causing interband transitions. 

In the description of PP experiments, we are interested only in the po- 
larization component propagating along the probe direction k^. For a weak 
probe, the nonlinear signal then arises from the linear response of the system 
in the presence of the pump, described by the time-dependent Hamiltonian 
H + Hp(t), to the probe-induced perturbation Hs{t). Within x^^\ the above 
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is true even for comparable pump and probe amplitudes. However, since, in 
contrast to H, the Hamiltonian H + Hp{t) does not conserve the number of 
carriers in each band separately, the calculation of the linear response func- 
tion is not practical. Therefore, we seek to replace H + Hp{t) by an effective 
Hamiltonian H{t) that does conserve the number of e-h pairs in each of its 
Hilbert subspaces (i.e., is "block-diagonal"). As derived in Ref. ||154|] , this 
can be accomplished in any given order in the pump field £p{t) and for any 
pulse duration by using a time-dependent Schrieffer- Wolff /Van Vleck canon- 
ical transformation ||158| , |160|| . Here it is sufficient to "block-diagonalize" the 
Hamiltonian H + Hp{t) up to the second order in Sp{t). The transformation 
that achieves this has the form e~'^'^e~^^[H + Hp(t)]e'^^e^'^ , where the anti- 
Hermitian operators Ti{t) and T2{t) create/annihilate one and two e-h pairs, 
respectively. 

We proceed with the first step and eliminate the single-pair pump-induced 
transitions in the time-dependent Schrodinger equation of the pump/bare- 
semiconductor system. 



d 



H-Hp{t) 



(100) 



This is achieved by substituting |\E'(t)) = C"^^*-*^ |x(t)) and acting with the 
operator e"^^*^*-* on the Ihs of Eq. (|100|) , 



.d_ 

''di 



H 



3^^^*^lx(t)) = e-^^W [/7,(t)]e^^W|x(t)). 



(101) 



(102) 



The anti-Hermitian operator Ti{t) has a decomposition 

fi(t) = 'p{t)e-'^^-'' - v\ty^^\ 

where V\t) and V{t) create and annihilate single e-h pairs, respectively. 
The effective Hamiltonian can be found from the condition that the terms 
describing single-pair interband transitions cancel each other in Eq. (|101|) . In 
Ref. [154], it was shown that the multiple commutators of V{t) with its time 



derivatives can be eliminated from Eq. (|101|) to all orders. By expanding 
Eq. ( |101|) and neglecting third or higher order terms in P, we obtain the 
following equation: 



' dt 



(103) 
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with initial condition oo) = 0. The formal solution is 



J —oo 



(104) 



Note that, since the Hamiltonian H conserves the number of e-h pairs and 
the optical transition operator f/^ creates a single e-h pair, V\t) also creates 
a single e-h pair. Furthermore, since both H and W conserve momentum, 
so does 'P'^{t). Eq. (|101[ ) then takes the form 



lx(()> 



H.c. 



where 



H{t)=H + ^ (fp(t)[p(t),f/t]+H.c.) 



\x{t)). (105) 



(106) 



is the sought time-dependent effective Hamiltonian that conserves the num- 
ber of e-h pairs and 'P^(t) is given by Eq. ( |103|) . Note that, since 'P^(t) is 



linear in the pump field the pump-induced term in H{t) [second term in 
Eq. ( |106| )] is quadratic (oc SpS*). The rhs of Eq. ( |105| ) describes the pump- 
induced two-pair transitions. These can be eliminated as well by performing 
a second canonical transformation, Ixi^)) = e^^^^^\^{t)). Following the same 
procedure, we use the anti-Hermiticity of T2(t) to decompose it as 



T2(t) = P2(t)e 



-2jkr,r 



2ikB-r 



(107) 



where Vlit) and 'P2{t) create and annihilate two e-h pairs, respectively. Sub- 
stituting the above expression for \x{t)) into Eq. ( |105| ) and requiring that all 
two-pair transitions cancel out, we obtain the following equation for Vlit), 



Mit) 



dt 



(108) 



Note that Vlit) only affects the PP polarization via higher order (S^) correc- 
tions, which are neglected here. However, it does determine the four-wave- 
mixing (FWM) polarization (see below). 

To obtain the condition of validity of this approach, it is useful to write 
down a formal solution ( p.04| ) of Eq. ( |103| ) in the basis of the N-hole many- 
body eigenstates, |aA^), with energies EaN, of the Hamiltonian H. Here the 
index a labels all the other quantum numbers, so that A^=0 corresponds to 
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the semiconductor ground state |0), denotes the one-pair states (exciton 
eigenstates in the undoped case, with a labehng both bound and scattering 
states), |a2) denotes the two-pair (biexciton in the undoped case) eigen- 
states, etc. In this basis, the solution of Eq. (|103|) can be written as 



l^^N ^\\V\t)\aN) 
(/?A^+ l|[/t|aA^) 



J — oo 



)p-r(t-t') 



(109) 



where we separated out the detuning Q and denoted tS.Ef 



al3 — J-^l3N+l — EaN- 

It can be seen that for resonant excitation (small fl) the rhs of Eq. ( |109|) 
is of the order of fiSptp. Thus, for short pulses, this parameter justifies the 
expansion in terms of the optical fields. For off-resonant excitation, this 
expansion is valid for any pulse duration provided that fiSp/Q < 1. Similar 
conditions can be obtained for the two-pair transition described by 7^2 • 

The nonlinear polarization Eq. (^) can now be expressed in terms of the 
linear response to the probe field: 



P(t) = /ie-*"-*(^(t)|f/|^(t)) = /ie-*"-'(<l>(t)|t/T(t)|$(t)), 
where, in the first order in £s(t), the state |$(t)) is given by 

|$(t)) = |$o(t)) - f dt'}C{t,t') \Ss{t')e^^^''+''^^^Ul{t') + H.c. 



Here }C{t,t') is the time-evolution operator satisfying 

z^/C(t,0=^(t)/C(t,t'), 



(110) 

l<^'o(t'))- 
(111) 

(112) 



Q-^zy/g ^i(*)^tgT'i(i)gT2(t) ig ^]2e (transformed) optical transition 



and 4(t) = e-^^Wg-t, 
operator. In Eq. (|111|) , |$o(^)) = /C(t, — oo)|0) is the time-evolved ground 
state |0); since H{t) conserves the number of e-h pairs, |$o(^)) contains no 
interband e-h pairs (in undoped semiconductors, it coincides with the ground 
state, |$o(^)) = |0)). /,From Eqs. ( |111| ) and ( |110| ), the polarization P(t) takes 
the form 



Pit) 



dt' ($0 (t) I Ut (t) lC{t,t') [Ss {t'Y^^ ■''+''^-^1 {t') 
~j L 

+ £:{t')e-'''^-'~'''-^UTit')]\%if)) 
- ( $0 (t') I [£s {t')e^' •'■+*"-"f/^ (t') 



+ S:{t')e 



'^^UT{t')]ic(t\t)UT{moit)) 



;ii3) 



73 



The above expression for the total polarization contains contributions prop- 
agating in various directions. To obtain the polarization propagating in a 
specific direction, one has to expand the effective-transition operator UT{t) 
in terms of Ti and T2. Using Eqs. ( |102|) and ( |107| ) and keeping only terms 
contributing to PP and FWM polarizations, we obtain [|154| 



Um = Ul(t) + Ul(t)e'^-'' + UFWM{t)e 
where (to lowest order in the pump field) 



-2i]s.r,r 



+ ■ 



(114) 



ulit) 

Ulit) 



u 



t 

FWM 



it) 



u^ + -[v{t),[u\v\t) 



+ 



v\t),\u\v{t) 



P{t),U^ 

u,vHt)],v^{t)]-\u,vl{t) 



(115) 



Here operators Ui{t) = U{t) and UpwMit) create one e-h pair, while Ulit) 
creates two e-h pairs (note that UpwM in Eq. ( |115| ) annihilates an e-h pair). 



Pump— probe polarization 



In order to extract the PP polarization from Eq. (|113D, one should retain 



only terms that are proportional to e*''" ''. Substituting Eqs. ( p.l4| ) into Eq. 
([TT3|) , we obtain Pk,(i) = e^^s-r-iuAt~r)[pm(^t) + P^'^\t)], where 



P«(t) = tfx' / dt'£,it'){Mt)\Umit,t')U\t')\%{t')), (116) 
J —00 

and 

p(2)(t) = z/x^ /* dt's,{t'){Mm2m{t,t')uKt')\Mt'))- (117) 

J —00 

Note that the above formulae apply to both undoped and doped semicon- 
ductors. 

Eqs. ( |116| )- ([TT7|) express the nonlinear PP polarization in terms of the 
linear response to the probe field of a system described by the time-dependent 
effective Hamiltonian ( 106 ). Such a form for the nonlinear response allows 
one to distinguish between two physically distinct contributions to the optical 
nonlinearities. Assuming that a short probe pulse arrives at t = r, consider 
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the first term, Eq. ( 116 ), which gives the single-pair (exciton for undoped 
case) contribution to the PP polarization. At negative time delays, r < 0, 
the probe excites an e-/ipair, described by state W {t)\^o{t)) ^ f/^|0); since 
the probe arrives before the pump, the effective transition operator coincides 
with the "bare" one [see Eqs. (|115|) and ( |104| )]. The first contribution to the 
optical nonlinearities comes from the effective Hamiltonian, H{t), governing 
the propagation of that interacting e-h pair in the interval (r, t) via the 
time-evolution operator /C(t,r). Note that since the pump pulse arrives at 
t = 0, for |r| 3> r~^, the negative time-delay signal vanishes. At t > 0, 
the e-h pair (exciton in the undoped case) "feels" the effect of the pump 
via mainly the transient bandgap shift, leading, e.g., to ac-Stark effect, and 
the change in the band dispersions (increase in effective mass/density of 
states), leading to enhanced e-h scattering (exciton binding energy in the 
undoped case). Note that H{t) also contains a contribution coming from the 
interactions between probe- and pump-excited e-h pairs, which are however 
perturbative in the doped case for short pulses or ofT-resonant excitation and 
lead to subdominant corrections. Importantly, the response of the system 
to the optically-induced corrections in H{t) takes into account all orders 
in the pump field, which is necessary for the adequate description, e.g., of 
the ac-Stark effect and the pump-induced changes in the e-h correlations. 
Indeed, although the pump-induced term in Eq. ( [10(j| ) is quadratic in £p, the 
time-evolution of the interacting e-h pair is described without expanding 
]C{t,T) in the pump field. On the other hand, the third-order polarization 
(x*-^-*) can be obtained by expanding /C(t, r) to the lowest order. The second 
contribution to the optical nonlinearities comes from the matrix element 
of the final state, {^olt)\U{t) in Eq. (|rT|). The latter, given by Eq. ([TT5D , 
contains the lowest order (quadratic) pump-induced terms which describe the 
Pauli blocking, pair-pair, and pair-FS interaction effects (exciton-exciton 
interactions in the undoped case |[154|| ). Note that the matrix element of 
the initial state contributes for positive time delays, i.e., if the probe arrives 
after the pump pulse. In this case, however, the pump-induced term in the 
effective Hamiltonian ( p.O(j| ) vanishes (since it lasts only for the duration of the 
pulse) so that for positive t > tp the PP signal is determined by the matrix 
elements rather than by H{t). If the probe arrives during the interaction of 
the system with the pump pulse (r ~ tp), both the effective Hamiltonian and 
the matrix elements contribute to the polarization. In this case, there is also 
a biexcitonic contribution [given by Eq. ( [117| )] coming from a simultaneous 
excitation of two e-h pairs by the pump and the probe. However, such a 
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biexciton state does not contribute to negative (r < 0) time delays. As 
can be seen from the above discussion, our theory separates out a number of 
contributions that play a different role for different time delays and excitation 
conditions. 

FWM polarization 

By extracting from Eq. ( |113| ) all the terms propagating in the the FWM 
direction, 2kp — k^, we obtain [|154|| (for delta-function probe Ss(t) = SsS{t — 
r)) 



2 c* 

s 



X 



($o(t)|f//C(t,r)4^^,(r)|$o(r)) - (t ^ r) 



;ii8) 



where the FWM transition operator Up^yn^{t) is given by Eq. ( |115D . It is con- 
venient to express Up^yj^it) in terms of the "irreducible" two-pair operator 
W^{t) = |pt2 _ pt^ satisfying 



.dW^{t) 
dt 



H{t),W\t) \ + ii£p{t)U^P{t). 



;ii9) 



In terms of W\t), the state f/pH/MWI^ol^)) in Eq. ([TT8|) can be presented 
as a sum of two- and one-pair contributions: 

uiwMit)mt)) = uw\t)\%{t)) - PmPmoit))- im 

For undoped semiconductors, |$o(^)) represents the ground state |0) of H. 
The operator W^^(t), being quadratic in the pump field (oc Sp), describes the 
simultaneous excitation of two interacting e-h pairs by the pump pulse and 
includes the biexciton and exciton-exciton scattering effects. Introducing the 
amplidutes Xk and $k as 



and 



^VMW|0) = Exk(t)al6Lk|0), 

k 



K:{t,t')ul^^,{t')\o) = E$k(t,t')4&lklo), 



;i2i) 



;i22) 
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where ^k{t, t') satisfies 

idt^^it.t') = ^(0|6_kak^(t)ai6lq|0)$q(t,O, (123) 
q 

with initial condition $k(^,^) = Xk.{^)i finally obtains for the FWM 
polarization 

PpWMit) = ie^(2kp-k.).r-^c.,(t+r)^(^ _ ^^^2^* t) - *k(T, t)]. (124) 

k 

Note that the third order polarization [corresponding to x*'^'*] is obtained by 
replacing H with H. 

B APPENDIX B 

In this appendix we clarify our convention for the time delay r and re- 
late it to the most commonly used conventions in PP and FWM. In the 

generic experimental configuration two laser pulses i?i(t)e*''^''"~*'^*^*~*^\ and 
£'2(^)6*'^^ '""*'^*'* *^'') respectively centered at time t = ti and t = t2 are inci- 
dent on a sample. Let us define At as, 

At^ti- 12 (125) 

and consider a FWM experiment where the signal is measured in the direction 
2k2 — ki. Then for a two- level-atom, the signal vanishes for At < 0, while 
for At > its amplitude, which decays with time as e~*/"^^, is determined 
by the Pauli blocking. In a system with Coulomb interactions (such as a 
semiconductor) a FWM signal is observed both for At < and At > 0. The 
At < signal is entirely due to the Coulomb interaction. 

In PP experiments, one usually chooses one pulse (the "pump") to have 
an amplitude Ep much larger than that of the other pulse (the "probe"), 
Es- As discussed in the text, a weak probe measures the linear response of 
the system (bare or dressed by the pump). If wc choose that Ep = E2, the 
pump induces coherent and incoherent populations when it arrives in the 
sample before the probe i.e. for t2 < ti. This is usually defined as "positive" 
time delay r = t2 — ti > in the PP literature. Note that r = —At, i.e., 
the "regular" sequence in PP experiments is the reverse of that of FWM 
experiments. For r < 0, the origin of the PP signal is that the probe creates 
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a linear polarization in the sample which lasts for time ~ ^2 = F"^ and, 
consequently, is scattered by polarization excited by the pump field. The 
signal observed for r < is therefore due to coherent effects. 

In FWM experiments, there is no restriction on the magnitude of the two 
incident fields E2 and Ei, which are often chosen to have amplitudes of the 
same order. Note however that, at the x*-^^ level, the FWM and pump-probe 
polarizations are linear in the Ei{t) field and thus the above linear response 
calculation applies even for comparable pump and probe amplitudes. 



C APPENDIX C 



In this Appendix we present the explicit expressions for the renormalized 
transition matrix elements in the presence of the pump excitation. The 
direct transition matrix element is given by 



k'<kp 



2 E KH{p',t) 

p'>kp 



n\(p'p;p';^)-n\(p'p;p';^) 



K/.(p'p;p;^) + ^e\(p'p;p;^) 



^ p'>kF 



X 



P|,(pp';p;i)+P,\(pp';p;i) 



(126) 



The first term on the rhs of the above equation describes the phase space 
filling contribution, while the rest of the terms are due to the mean field 
pair-pair and pair-FS interactions. 

The pump-induced indirect transition matrix element is given by 



Veh{Kt)-Veh{p + p' -K t) 
Veh{p,t)+Veh{p',t) 
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X 



n\(pp'; p; t) + nt(pp'; p + p' - k; t) 



+ 7'e;.(p + p'-k,t) 



P|,(k,p + p'-k;p';i) 



7',\(k,p + p'-k;p;i) 



P,\(k, p + p' - k; p'; + P|,(p, p'; k; 



- Kf,{p, p'; p + p' - k; - 7',\(k, p + p' - k; p; t) 

-Veh{p',t)V2{Kp + p'-k;p;t). 
The effective e-h potential is given by 



Vehiq.] kk'; t) = v{q) - ^ £p{t) 



V!d^ + q,k';k' + q;t) 

+ V^h{k,k' + q;k';t) 
+ 7'-(k,k' + q;k';i) 

+ P,\*(k + q,k';k' + q;t) 



and tfie effective e-e potential is given by 



Vee{q;kk';t)^v{q) + f^Sp{t) 



P|,(k + q,k'-q;k';i) 

-7'|Jk + q,k'-q;k;i) 
+ 7'|,*(k,k';k'-q;t) 

-P,t(k,k';k + q;t) 



(127) 



(128) 



(129) 
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Figure 1: The effect of the pump pulse on the band dispersion, £\i{t) = 
£k(i) + £^-k(^)) of the "pump-dressed" system. Sohd Une: "bare" dispersion 
{t = — oo). Dashed hne: pump-renormahzed dispersion (t — 0). The bands 
become "heavier" for the duration of the pump. 
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Figure 2: (a) The e-h scattering processes that contribute to the rhs of Eq. 
(|39|). Full lines correspond to s(p, k, t) and thin lines to V. The diagrams 
describe (from top to bottom) (i) Born scattering of a FS electron, (ii) FS 
electron ladder diagrams, (iii) FS hole ladder diagrams, and (iv) Nonlinear 
vertex corrections due to the dynamical FS response, (b) Scattering processes 
state that determine the time- and momentum-dependence of the effective 
e-h potential V{p, t) and lead to the unbinding of the HFA bound state. 
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Figure 3: Linear absorption resonance lineshape for FES (solid curve) com- 
pared to the HFA (dashed curve), calculated with g — 0.4 and F = O.IE'f. 
The HFA resonance position was shifted for better visibility. 
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Figure 4: The absorption spectra for the FES (a) compared to the HFA 
(b). Sohd curves: hnear absorption. Dashed curves: nonhnear absorption. 
The curves were calculated with g = 0.4 and F = O.lEp at time delay 
T = — 0.ir~^ = ~^p/2 and pulse duration tp — 2.0Ep^. The nonlinear 
absorption lineshapes exhibit bleaching, resonance blueshift, and gain below 
the absorption onset that differ in the two cases. 
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Figure 5: Differential transmission linesliape for various time delays cal- 
culated with g = 0.4, T = O.lEp, and tp = 2.0Ep^. (a) Long time delay, 
T = — 1.5r~^ — —Ib.QEp^ . For the FES, the oscillations in the differen- 
tial transmission spectra are reduced (solid curve) as compared to the HFA 
(dashed curve), (b) Short time delay, r = -O.lTs = -tp/2 = -1.0^;^^ For 
the FES, the differential transmission spectrum is asymmetric (solid curve) 
as compared to the symmetric lineshape for the HFA (dashed curve). The 
above curves were shifted for better visibility. 
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Figure 6: (a) Nonlinear absorption resonance bleaching evaluated at the 
instantaneous peak frequency as function of time delay for the FES (solid 
curve) compared to the HFA (dashed curve). The curves were calculated with 
g = 0.4, F = O.lEp, and tp = 2.0Ep^. The time-dependence of the pump 
pulse is also presented for comparison (dotted curve), (b) Same calculated 
using the rigid band shift model, Eq. (^). 
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Figure 7: The enhancement of the nonhnear absorption resonance of the 
FES (dashed curve) vs. hnear absorption (sohd curve) at long time delay 
T = — 1.5r~^. The curves were calculated with g — 0.4, F = Q.lEp-, and 
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Figure 8: Nonlinear absorption resonance blueshift as function of time delay 
for the FES (solid curve) compared to the HFA (dashed curve). The blueshift 
is significantly weaker for the HFA. The curves were calculated with g = 0.4, 
r = O.lEp, and tp — 2.0Ep^. The time-dependence of the pump pulse is 
also presented for comparison (dotted curve). 
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Figure 9: Resonance bleaching as function of time delay for the FES (a) 

compared to the HFA (b) for different strengths of the e-h interaction: g=0.4 
(solid curve) and g=0.3 (dashed curve). The curves were calculated with 
r = O.lEp, and tp = 2.0Ep^. The time-dependence of the pump pulse is 
also presented for comparison (dotted curve). 
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Figure 10: The effect of tlie pump-induced renormalization of tlie band dis- 
persions on tlie e-h interactions. Tlie function F{uj,t), given by Eq. (|60D, 
for the FES (a) compared to the HFA (b) in the presence (dashed curve) and 
absence (solid curve) of the pump pulse. The curves were calculated with 
g = 0.4, r = O.lEp, and tp = 2.0Ep\ 
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Figure 11: The effect of the pump-induced renormaUzation of the band 
dispersions on the resonance strength. The nonhnear absorption spectrum 
for the FES (a) compared to the HFA (b). Sohd curves: Linear absorption. 
Dashed curves: Nonhnear absorption. Dotted curves: Nonhnear absorption 
for a rigid band shift only. The curves were calculated with g — 0.4, F = 
O.lEp, and tp = 2.0Ef^. 
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Figure 12: (a) Calculated PP spectra with (a) A = l.TEp, (b) A = l.QEp, 
and (c) A = l.&Ep (HFA), for short pump duration tpEp/% ~ 2.0, and 
negative time delays tT /U — —2.0 (lowest curve), —1.2, —0.6, —0.4, —0.2, 0, 
and hnear absorption spectrum (upper curve). Inset: schematic plot of the 
energy spectrum of the two-subband QW (right) and absorption spectrum 
(left). 
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Figure 13: Calculated absorption spectra at positive time delays for 
nanoparticles with (a) R — 5 nm, (b) R — 2.5 nm, and (c) R — 1.2 nm. 
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Figure 14: Calculated differential transmission spectra at positive time 
delays for nanoparticles with (a) = 5 nm, (b) R = 2.5 nm, and (c) R = 1.2 
nm. 
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Figure 15: Temporal evolution of the differential transmission at frequencies 
close the SP resonance for nanoparticles with (a) R — 5 nm, (b) R — 2.5 nm, 
and (c) R — 1.2 nm. (d) Measured time-resolved pump-probe signal. 
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